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ABSTRACT 

We study how optical galaxy morphology depends on mass and star formation rate 
(SFR) in the Illustris Simulation. To do so, we measure automated galaxy structures 
in 10808 simulated galaxies at z = 0 with stellar masses 10^-^ < M^/Mq < 10^^-^. We 
add observational realism to idealized synthetic images and measure non-parametric 
statistics in rest-frame optical and near-IR images from four directions. We find that 
Illustris creates a morphologically diverse galaxy population, occupying the observed 
bulge strength locus and reproducing median morphology trends versus stellar mass, 
SFR, and compactness. Morphology correlates realistically with rotation, following 
classification schemes put forth by kinematic surveys. Type fractions as a function 
of environment agree roughly with data. These results imply that connections among 
mass, star formation, and galaxy structure arise naturally from models matching global 
star formation and halo occupation functions when simulated with accurate methods. 
This raises a question of how to construct experiments on galaxy surveys to better 
distinguish between models. We predict that at fixed halo mass near disc- 

dominated galaxies have higher stellar mass than bulge-dominated ones, a possible 
consequence of the Illustris feedback model. While Illustris galaxies at r\j lO^^Mo 
have a reasonable size distribution, those at ^ have half-light radii larger 

than observed by a factor of two. Furthermore, at ^ 10^^‘^-IO^^Mq, a relevant 
fraction of Illustris galaxies have distinct “ring-like” features, such that the bright 
pixels have an unusually wide spatial extent. 

Key words: galaxies: structure — galaxies: statistics — galaxies: formation — meth¬ 
ods: numerical 


1 INTRODUCTION 

Cosmological simulations of galaxy formation are poised to 
address long-standing challenges in extragalactic astronomy. 
They may achieve this not only by the ordinary means of 
testing theories of galaxy astrophysics with numerical exper¬ 
iments, but also by enabling the direct modeling of observa¬ 
tions. In observed distant galaxies, often our only dynami- 


* Hubble Fellow 


cal information is a snapshot of the stellar orbits traced by 
light, called the galaxy’s morphology. Despite copious data, 
the origins of the distribution of galaxy morphologies remain 
debated. The presence of a bulge-like light prohle correlates 


tightly with a lack of star formation (e.g., Kauffmann et al 
2^[Bell etal.|2012 |Woo et al.|2015| ), but we know neither 


the complete physics of quenching nor the assembly histo¬ 
ries of bulges as a function of halo mass or environment. 
Although rare, galaxy mergers can play an important role 
in massive galaxies (e.g.. Bell et aL]|2006 Lotz et al.|2008b 
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2011). However, their rates and consequences are largely un¬ 


known at z > 1. 

Analyses of galaxy morphologies typically span two 
broad properties: galaxy structure and disturbed appear¬ 
ance. Important methods for conducting these analyses are 
visual inspection and automated algorithms, both of which 
have matured as the quantity of high resolution galaxy data 
increases. However, the results of such analyses are for the 
most part empirically motivated: they determine whether a 
galaxy appears to be bulge- or disc-dominated and whether 
it appears to be experiencing a merger, interaction, or dis¬ 
turbance. We often define these rigidly, using wide bins that 
may be ineffective for learning how galaxies formed. 

In principle, simulations can link observed morphology 
with underlying physical processes in powerful ways. The 
principal idea is to apply stellar population synthesis tech- 


niques (e.g.. 

Tinsley|I968j Tinsley V Gunn3l976 Tinsley V 

Larson|I978 

Gunn et al.|I98I) to galaxy evolution models. 


For example, analytic and semi-analytic models make uni¬ 
fied predictions for observed fluxes, dynamics, and structure 


based on simple physical prescriptions (e.g., Henriques et al. 


2015 Lu et al. 2014). Indeed, these calculations now form the 


basis for “mock observatories” in which we predict galaxy 
properties, including morphology, in realistic synthetic imag¬ 
ing surveys and release these to researchers through ad¬ 


vanced databases (e.g., Overzier et al. 2013 Bernyk et al. 
2014). However, by design they do not trace the full dynam¬ 


ics of matter in galaxies, especially rare and subtle features 
such as tidal tails, disturbed morphologies, and rapid mor¬ 
phological transformations. In contrast, using hydrodynam- 
ical simulations with sufficiently high spatial resolution and 
enough physics to predict chemically enriched stars, gas, and 
dust, we can assign spectral energy distributions (SEDs) to 
their components and predict how they would be observed in 
imaging surveys. Moreover, because such simulations model 
directly the gravity of baryons and dark matter, they also 
capture processes such as local collapse and angular momen¬ 
tum transport which are essential for setting galaxy mor¬ 
phologies. 

Until recently, the direct simulation approach has been 
limited either to very small samples of cosmologically as¬ 
sembling galaxies or to hand-crafted models of individual 
galaxies or groups. However, significant progress has been 
made toward achieving direct, statistically relevant predic¬ 
tions of galaxy formation. By including significant feedback 
from supernovae, massive star formation, and supermassive 
black holes, the realism of resulting galaxy models has con¬ 
tinued to improve (e.g., Crain et al.|20Q9 Guedes et al. 2011 


Ceverino et al.|2014 Cen|2Q14 ). These advances owe in part 

to preceding generations of semi-analytic modeling which 
outlined the broad requirements for the impacts of feedback 


on massive galaxies (e.g., Croton et al. 2006 Somerville et al. 
2008). Moreover, refinements in numerical methods have en¬ 


abled the field to conduct very large hydrodynamical simula¬ 
tions containing tens of thousands of galaxies whose internal 


dynamics are at least partially resolved (e.g.. 

jSchaye et al. 

20I01JKhandai et al.||20I5 

IDubois et al.||20I4| 

IVogelsberger 

et al.|20I4b| iSchaye et al. 

2014). Since they predict the dis- 


tribution and appearance of rare events, and directly link 
the physics of star formation quenching with the evolution 
of galaxy structure, these simulations increase the informa¬ 
tion learned from galaxy surveys. We can, in principle, dis¬ 


card empirically motivated and coarsely binned classification 
schemes in favor of explicit ones. 

As a first step, in this paper we study the morphology 
of 10808 galaxies at z = 0 from the Illustris Project. We 
measure each galaxy from four viewing directions in four fil¬ 
ters, for a total of ~ 1.7 x 10^ synthetic images. By doing 
so, we seek to enable direct comparison between simulated 
and observed galaxy populations and to prepare for future 
studies exploiting these models to interpret survey data. In 
Sectionwe describe our methods for creating realistic syn¬ 
thetic observations from ideal images. In Section we de¬ 
termine how these galaxy images compare with observations 
using non-parametric image diagnostics, and show how Il¬ 
lustris galaxy models occupy the space of stellar mass, SFR, 
and morphology. Sectionj^links optical morphology with the 
kinematic structure of the simulated galaxies. In Section 
we explore correlations of morphology with galaxy size, the 
mass of supermassive black holes, and halo mass. We will 
make our non-parametric morphology measurements avail¬ 
able as online supplementary tables (see Appendix E . We 
have released our synthetic image^ and publicly releasing 
our code for adding realism to model image^ 

For all calculations, we adopt the cosmology used for the 
Illustris simulation, a ACDM cosmology with Qrn = 0.2726, 
Qa = 0.7274, Qb = 0.0456, erg = 0.809, n, = 0.963, and 
Ho = lOO/i kms“^ with h = 0.704, consistent with the 
WMAP-9 measurements ( Hinshaw et al.||2013 ). 


2 MOCK IMAGE ANALYSIS 

In this section we briefly outline our methods to create realis¬ 
tic mock images and measure morphologies in observational 
units. 


2.1 Simulations and Ideal Images 


The Illustris Project consists of hydrodynamical simulations 
of galaxy formation in a periodic volume of size (106.5 Mpe)^ 
( Vogelsberger et al.|20I4a]b Genel et al.|20I4 ). This volume 
was simulated with and without baryons, and at several res¬ 
olution levels; unless otherwise noted, here we present results 
from the simulation with the highest resolution (2 x 1820^), 
which resolves baryonic matter with mass 1.26 x 10^ Mq and 
employs a plummer equivalent gravitational softening length 
of 710 pc at z = 0. Using the Arepo code for gravity and gas 
dynamics ( Springel|20I0 ), these simulations evolve simulta¬ 
neously the gas, stars, and dark matter components from 


cosmological initial conditions (Vogelsberger et al. 2012 


Keres et al. 2012 Sijacki et al. 2012). The Illustris simulation 


galaxy physics models consist of primordial and metal line 
cooling, star formation, gas recycling, metal enrichment, su¬ 
permassive black hole (SMBH) growth, and gas heating by 
feedback from supernovae and SMBHs (Vogelsberger et al. 


2013). Parameters were chosen to match the z = 0 stellar 


mass and halo occupation functions, plus the cosmic his¬ 


tory of SFR density. Torrey et al. (2014) showed that this 


1 www.illustris-project.org 
^ https://bitbucket.org/ptorrey/sunpy 
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parameterization reproduces the observed stellar mass func¬ 
tion, SFR-Mass main sequence, and the Tully-Fisher rela¬ 
tion, from z = 3 to z = 0. Sijacki et ah (2015) described 
the evolution of SMBHs in Illustris and their co-evolution 
with the galaxy population. The Illustris simulations are 
publicly available at www.illustris-project.org/data ( [Nelson 
et al.||2015b ) 


We generate noiseless high-resolution images with an 


built around the Sunrise 


et al. 


2010 



code ([Jonsson 

2006 

Jonsson 


Jonsson & Primack 2010). T15 published ide¬ 


alized images to enable studies focused on data from dif¬ 
ferent observatories with the same simulated dataset. For 
example, Wellons et al. (2015) used such tools to demon¬ 
strate how high-redshift compact massive galaxies would be 
observed with the Hubble Space Telescope. Similarly, jMoodyj 


et al. (2014a) created and analyzed mock HST images from 


very high-resolution cosmological simulations of clumpy disc 
galaxies ( Ceverino et al.|20T4 ). We describe our methods for 
adding realism and our example model observatory settings 
tuned for low-redshift Illustris galaxies in Section [2.2| 

To create the images, we begin by segmenting the sim¬ 
ulation output into separate hies, one for each galaxy (sub¬ 


halo) identihed with the SUBFIND code (Springel et al. 


2001 ) 


et al. 


, which uses a friends-of-friends algorithm (e.g., Davis 


1985) with a linking length of 0.2 times the mean parti¬ 


cle separation to identify bound dark matter halos, and then 
applies a secondary linking stage to associate baryonic mass 
to these halos ( Dolag et al.|2009' ). In this paper, we analyze 
only the mass associated with individual subhalos, ignoring 
nearby subhalos such as other members of the parent halo. 
This allows us to separate galaxies into images and cleanly 
measure structural parameters. However, with this approach 
we are unable to directly study signatures of galaxy mergers 
such as multiple nuclei, since the halo finder may separate 
objects that would be combined in real image-based segmen¬ 
tation maps. Moreover, the procedure used here occasionally 
leaves spurious faint surface brightness features from satel¬ 
lite galaxies; these do not affect the global optical morphol¬ 
ogy of the bright pixels, but could affect diagnostics of faint 
or subtle stellar halo features. For studying galaxy merger 


signatures with simulations like Illustris (e.g., Rodriguez- 
Gomez et al.|2015 ), we will use images containing all nearby 


mass. 

With Sunrise, we assign to each star particle a spec¬ 
tral energy distribution (SED) based on its mass, age, and 
metallicity using stellar population models by |Bruzual~fcl 
Chariot (2003) with a jChabrie^ (2003) initial mass func¬ 


tion. Each star particle radiates from a region with a size 
directly proportional to the radius that encloses its 16 near¬ 
est neighbors: the emission region is smaller (larger) in areas 
of high (low) stellar density. This adaptive smoothing tech¬ 
nique leads to much more realistic surface brightness maps 
in regions of low stellar density; see T15 for full details. 
This step preserves the light distribution in bright regions 
of well resolved galaxies, and prior to any analysis, we con¬ 
volve with a point-spread function (PSE) larger than the 
typical spacing between star particles near the centers of 


^ Sunrise is freely available at 

https://bitbucket.org/lutorm/sunrise 
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galaxies. Therefore, the results of our study do not depend 
strongly on this assignment. However, certain measures of 
dumpiness, such as Gini alone, are sensitive to individual 
bright star particles and therefore care must be taken when 
interpreting the original high-resolution images regardless of 
the light assignment method. In Appendix A, we used a rep¬ 
resentative subset of galaxies to compare our choice of star 
particle sizes with the most common alternative, a constant 
radius equal to stars’ gravitational softening length, and hnd 
that the morphology diagnostics measured here are identical 
to within 1%. 


We then project these SEDs with Sunrise without per¬ 
forming dust radiative transfer (RT) through the ISM. We 
skip the RT stage in order to avoid over-modeling galax¬ 
ies with ISM resolution insufficient to gain substantial ac¬ 
curacy by performing the RT , and because the primary 
focus of the present work is on the morphologies of mas¬ 
sive > 10^'^Mq) galaxies at z = 0, which we do not 
believe to be especially senstive to dust attenuation, on av¬ 
erage. Therefore, in this study we have used Sunrise only 
as a tool to assign stellar population models to star parti¬ 
cles, project these quantities in several directions, and cre¬ 
ate synthetic images in arbitrary hlters. While this may be 
slightly more computationally intensive than required for 
the present work, we viewed the advantages of starting with 
a widely used open-source code as the backbone of our anal¬ 
ysis pipeline, with the option of using dust RT in future 
studies or in higher resolution simulations, to outweigh the 
disadvantages. 


To assist our analysis, we also project surface density 
maps of physical quantities onto the same pixels, using the 
options provided by Sunrise. These maps include SER, 
stellar mass, gas mass, stellar metallicity, gas metallicity, 
luminosity- and mass- weighted stellar age, and gas tempera¬ 
ture. In this paper, we neglect the important effect that dust 
attenuation will have on our synthetic images. However, we 
have performed extensive tests derived from these 2-D mass 
maps, from which we derived a resolved slab dust model by 
expanding on techniques to estimate dust on semi-analytic 
models (e.g., De Lucia Blaizot|2007 Kitzbichler & White 


2QQ7||Guo fc White|20Q9||Guo et al.|2QII||20I2| ). We found 


that the primary results and correlations presented in this 
paper are qualitatively unaffected by dust modeled in this 
way, so we have chosen to reserve a comprehensive explo¬ 
ration of this topic to a future paper. 


The unattenuated rays travel until they exit the grid or 
enter one of four viewing apertures (“cameras”). The out¬ 
put of this step is the SED at each of 512x512 pixels in 
each camera. Erom these data cubes. Sunrise creates raw 
mock images by integrating the (optionally redshifted) SED 
in each pixel over a set of common astronomical hlters, from 
the UV through IR. In this paper, we perform this filter syn¬ 
thesis in the rest frame of each galaxy. The spatial extent 
of each image is set to ten times the 3D stellar half-mass 
radius, and therefore the physical pixel scale varies. We use 
pixel sizes of roughly 100-300 pc, sufficiently small to sim¬ 
ulate SDSS, Pan-STARRS, and LSST images of sources at 
z > 0.02 regardless of the limited resolution of the simula¬ 
tion, which is a separate constraint that we consider below. 
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Figure 1. Example g-r-i images of 2 ; = 0 Illustris galaxies at each of three stages of our image realism procedure, arranged by increasing 
stellar mass from bottom (M* ~ 10^® Mq) to top (M* ~ 10^^ Mq) and selected based on their location in the G-M 20 plane (as in 
Section [ 3 ^ . For displaying, we pretend the galaxy was observed in SDSS at 2 ; = 0.05 with a pixel size of 0.5 arcsec. Left: Ideal output 
image directly from Sunrise. Middle: We convolved the ideal image for each filter by a Gaussian PSF with FWHM= 2.0 arcsec. Right: 
To the PSF-convolved images, we add random cutouts from SDSS DRIO ( |Ahn et al.||2014| >, downloaded from data.sdss3.org/mosaics. 
The remainder of this paper analyzes measurements of synthetic images that are most like those in the middle column, as described in 
Section 
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2.2 Realistic Images 

We convert the noise-free, ideal galaxy images (T15) into 
realistic synthetic images using the following procedure. 
First, we convolve each high-resolution image with a Gaus¬ 
sian point-spread function (PSF) with full-width at half¬ 
maximum (FWHM) of 1.0 kpc. Then, we re-bin the images 
to a constant pixel scale of 0.24 kpc, which is approximately 
1/3 of the Illustris stellar gravitational softening. These pa¬ 
rameters correspond roughly to 1 arcsec seeing for observa¬ 
tions of a source located at = 0.05, where the physical scale 
is roughly 1 kpc/arcsec (angular-size distance Da ~ 200 
Mpc, luminosity distance Dl ^ 220 Mpc). Thus our realis¬ 
tic images can be thought of as roughly appropriate compar¬ 
isons to many sources in the Sloan Digital Sky Survey main 


& White 2007) will become very useful as the volumes of 


galaxy sample (Strauss et al. 2002). They also correspond 
roughly to an HST WFC3 survey of sources at z ~ 0.5, 
where the scale ~ 6 kpc/arcsec implies that our PSF FWHM 
corresponds to ~ 0.17'k 

We save FIT^ images, either in the Sunrise output 
units {yV ! vn! vn? I Sr) ^ which is a surface brightness inde¬ 
pendent of distance, or in nanomaggies (1 nmy 3.631 x 
10“®Jy) assuming our fiducial setting z = 0.050. To sim¬ 
plify analysis, we also compute AB absolute magnitudes, 
AB absolute zero-points, original camera distances and pixel 
scales, and the new implied apparent magnitude at the cam¬ 
era distance. We store these important metadata in the im¬ 
age headers. Thus the conversion from the Sunrise output 
can be uniquely specified, and synthetic image fluxes can be 
recomputed for any assumed distance. From these FITS files, 
we then create files containing colour-composite images. 

Finally, we add sky shot noise such that the average 
signal-to-noise ratio of each galaxy pixel is 25. We assume 
this sky shot noise is a Gaussian random process indepen¬ 
dently applied to each pixel. Thus we are assuming that each 
model galaxy is strongly detected, eliminating biases from 
potentially noisy morphology measurements. 

For future visual classification projects, we also prepare 
images for classification by the Galaxy Zoo project (GZ, 
e.g., Lintott et al.|[2Q0^ in SDSS g, r, and i filters. Us¬ 
ing our initial radius (rp) measurements as defined in Sec- 
tion |2.3| we re-bin our SDSS-like FITS images to a new pixel 
scale (0.008 x Vp) and create images with a fixed pixel count 
(424 X 424). These choices are such that the galaxy extent de¬ 
fined by 2rp always subtends ~ 2/3 of the linear image size, 
enabling fair visual classifications as a complement to our 
fixed-scale non-parametric measurements below. For such 
visual classification projects, we also add real SDSS back¬ 
ground images to create fully synthetic ugriz galaxy images. 
To accomplish this, we first downloaded mosaics from the 
SDSS DRIO ( Ahn et M^|2014 ) Science Archive Server with 
the mosaic web tool (data.sdss3.org/mosaics). From these, 
we randomly select a region of an appropriate size for each 
synthetic image, assuming the galaxies are at z = 0.05, and 
add it to the simulated galaxy image. We demonstrate these 
steps in Figure This is a simplification from complete 
image simulations of self-consistent lightcones drawn from 
the simulation volume (e.g., Overzier et al.|20T3 Henriques 


et al.|2012 |. We have created several examples of these sim¬ 
ulated helds from Illustris, and such techniques (Kitzbichler 


such simulations grow. 


2.3 Structural Measurements 

We measure non-parametric morphologies by using code 
originally developed for idealized merger simulations and 
also applied to galaxy surveys (Lotz et al. 2008a|b 2011 
e.g.,). From each image, we characterize the light profile 
with a Petrosian radius rp, half -light radius R 1 / 2 , Goncen- 
tration (C), Gini (G), and M 20 ( Gonselice et al.| 2003 Lotz 
et al.||20Q4| , defined below. Our code also measures merger 
and disturbance indicators. Asymmetry and the newly pro¬ 
posed MID merger statistics ( Freeman et al.|[2013 ), but in 
this paper we focus on the above simple estimates of galaxy 
structure. We will return to diagnostics of mergers and dis¬ 
turbances in a future paper. 

We define the Petrosian radius or semi-major axis Vp 
such that the mean surface brightness in an elliptical annulus 
with semi-major axis Vp equals 0.2 times the mean surface 
brightness within this ellipse, following Lotz et al.| ( 2004| ). 
Here we also compute an elliptical half-light radius R 1/2 to 
characterize galaxy sizes, assuming that all of the galaxy’s 
light is contained within an ellipse with semi-major axis 1.5 x 


We compute the concentration parameter C (Bershady 
et al.|2000| l: 


rso 


C7 = 51ogio — 


r2o 


( 1 ) 


where rso and r 2 o are circular apertures containing 80% and 
20% of the total flux within the ellipse with semi-major axis 
1.5rp (Gonselice et al. 2003) of the galaxy center defined 


by minimizing the Asymmetry parameter (Abraham et al. 
19^. 


Gini’s coefficient, G, measures the inequality in flux 
value among a galaxy’s pixels, varying from 0 (all pixels 
equal flux) to 1 (one pixel contains all flux). First used by 
Abraham et al.| ( |2003| to characterize galaxy light profiles, 
G correlates with G but does not depend on the location of 
the brightest pixels. Hence it is sensitive not only to concen¬ 
trated spheroids but also to galaxies with multiple bright 
regions. For a discrete population, Glasser (1962) showed 
that G can be computed as: 


G = 


1 


\Ii\n{n - 1 ) 




( 2 ) 


where we have n pixels with rank-ordered absolute flux val¬ 
ues IA I, an d |A| = IAI /n, the mean absolute flux value. 
We follow Lotz et al.| ( |2Q04| in correcting G using abso¬ 
lute values to mitigate the effect of noise-induced negative 
fluxes. This procedure recovers the true G when S/N > 3 
per galaxy pixel, which is true by construction for all of the 
galaxy images we prepared in Section [2.2| 

M 20 measures the second-order spatial moment of a 
galaxy’s bright pixels contributing 20% of the total light, 
relative to its total moment ( Lotz et al.||2004 ): 


^ Flexible Image Transport System, fits.gsfc.nasa.gov 


M 20 = logio for V 7i < 0.27tot, 

Mtot 


( 3 ) 


© 0000 RAS, MNRAS 000 , 000-000 




























































6 G.F. Snyder et al. 


0.7 


0.6 

'E 

b 0.5 


0.4 


0.3 


Figure 2. Gini-M 2 o as a function of star formation as reflected hy U — B colour for galaxies at z ~ 0.3. Left: data from the Extended 
Groth Strip (EGS) survey compiled by |Lotz et al.| ( [2008b| >: roughly rest-frame B-band. Right: mock data from T15 at 2 : = 0.33: rest- 
frame r-band. Here we show simulated galaxies dX z = 0.33 to roughly match the volume-limited EGS sample, but all other flgures show 
simulated galaxies at 2 ; = 0 only. Red (blue) contours encircle 68% and 95% of the galaxies with U — B > 1 {U — B < 1). We And that 
simulated galaxies of a given colour have roughly the right optical shape, and Illustris produces a large population of early-type objects. 
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where 


I-0.693M20 + 4.95G - 3.85|, and 


n n 

Mtoi = '^Mi ='^Ii [(®i - + iVi - Vc)^] , (4) 

i i 


F(G,M2o) 


|F| G^0.14M2o +0.778 

-\F\ G < O.I 4 M 20 + 0.778 


and Xc, Vc are the 2-D spatial coordinates of the galaxy cen¬ 
ter, defined to minimize Mtot- For computing G and M 20 , 
we define a galaxy’s pixels following the segmentation pro¬ 


cedure by Lotz et al. (2004). 

In Figures and we divide the G-M 20 plane into 
three regions corresponding to early types, late types, and 
mergers, based on comparisons with low-redshift visual clas¬ 
sifications ( Lotz et al.|2004 ). These classifications are meant 
to be loose guidelines and are not to be strictly inter¬ 
preted. For future reference, we define a G-M 20 “bulge statis¬ 
tic”, which depends on an object’s location in this diagram 
and correlates with optical bulge strength. This will serve 
as a rough automated assessment of morphological type. 
Specifically, we define F as five times the point-line dis¬ 
tance from a galaxy’s morphology point to the line pass¬ 
ing through (Go, M 2 o,o) = (0.533, —1.75) and parallel to the 
Lotz et al. ( 2004| early-type/late-type separation line, which 
has a slope m = 0.14 in the space of (G, M 20 ). We chose the 
scaling factor 5 so that the resulting values occupy a con¬ 
venient range (—2 < F < 1). Starting from the point-line 
distance formula in two dimensions: 


d = 


|aM 2 o + 5G + c| 

(a2 + b2)l/2 


(5) 


where d is the distance from a point to the line G = 
— {a/b)M 2 o — (c/b). We let 5 = 1 and set a = 
—m = —0.14, allowing us to solve for c = —5(Go + 
ttM 2 o,o) = —0.778, which defines the desired line. We set 
the sign of F so that positive (negative) values indicate 
bulge-dominated (disc-dominated) galaxies. Thus, |T| = 


corresponding to the “G-M 20 bulge statistic” annotation to 
Figure]^ and shown in panel (d) of Figure]^ For most galax¬ 
ies, this diagnostic adds little new information beyond M 20 
(or G or Sersic index). However, F is less sensitive to dust, 
mergers, and other disturbances that move galaxies in a 
roughly perpendicular direction away from the main G-M 20 
locus. F traces quenched galaxies similarly well, if not a lit¬ 
tle better, than other structural parameters, and F is closely 
related to the dominant principal component of a suite of 
ten non-parametric structural diagnostics (M. Peth et al. in 
prep.). 

With these, we can automatically characterize the Il¬ 
lustris galaxy morphologies and compare them directly to 
measurements from real galaxies. In Appendixwe present 
the morphology catalog for Illustris at 2 : = 0. 


3 SIMULATED MORPHOLOGIES 

Hydrodynamical models have only recently begun to repro¬ 
duce basic statistics of the galaxy population, such as the 
galaxy stellar mass function and global star formation his¬ 
tory. Without these ingredients, it has been impractical to 
compare morphology surveys to simulations in interesting 
numbers. However, given recent achievements, it is worth¬ 
while to evaluate statistically the predicted morphologies 
and how they depend on mass and star formation. The ex¬ 
tent to which they match may provide clues to problems 
with the assumed galaxy formation physics, or may reveal 
the relative sensitivity of observed galaxy populations to the 
assumed physics. 
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In this section, we examine our morphology measure¬ 
ments from Illustris at ~ 0. This approach will also sup¬ 
port studies of how these morphologies change with time or 
would change if the galaxies were assumed to be at a greater 
distance, with shallower images, or with various observato¬ 
ries. 


For brevity, we have neglected dust. Information about 
the metal and gas content of the simulated galaxies can be 
used to predict the effects of dust, for example using the 
2 -dimensional projections of metal density described in Sec¬ 
tion We tested this technique, and found that the most 
important effect of dust is to lower the G value, and raise 
the M 20 value of simulated galaxies, especially where G is 
very large (for example, in the “mergers” region). In fact, 
the observed distribution (Figure]^ implies that some kind 
of attenuation is required to bring some high-G, high-M 2 o 
simulated galaxies into better agreement with observed ones. 
We have found that these general trends are indeed repro¬ 
duced by the simple dust treatments we have tested (follow¬ 
ing, e.g., De Lucia & Blaizot 2007). Specifically, in rest-frame 
u filter images of galaxies with G > 0.6, the effect of dust 
is to reduce G by ~ 0.05, on average. In other words, it 
tends to equalize the distribution of flux among the galax¬ 
ies’ pixels by obscuring the brightest star-forming regions. In 
galaxies with G ~ 0.5, M 20 ~ — 2 (typically relatively mas¬ 
sive composite disc+bulge systems), dust serves a similar 
role by increasing M 20 by ~ 0.2 on average, i.e. attenuat¬ 
ing the bright cores and causing the light prohle to appear 
slightly more extended. However, these effects largely do not 
change the coarse structural type that we would assign to 
the simulated images. Moreover, along the main locus of 
points, the average effect of dust imparts no systematic bias 
on the morphology measurements. These results are con¬ 
sistent with simulated image analyses by Lotz et al. (2008a) 
and Snyder et al. (2015). We will describe this dust modeling 
and related issues in a separate paper. 


3.1 Comparison with data 

In Figure we compare the 1200 most massive Illustris 
galaxies at 2 : = 0.33 with a volume-limited sample of galax¬ 
ies in HST imaging of the Extended Groth Strip (EGS) at 
0.2 < z < 0.4 ( Lotz et al.|2008b ). These samples have a sim¬ 
ilar H-band absolute magnitude upper limit {B < —19.5). 
Eor a direct comparison in Eigure we used simulated 
galaxies at z = 0.33, whereas the rest of this paper dis¬ 
cusses z = 0 only. We added realism to the simulated im¬ 
ages at z = 0.33 so that their resolution and noise levels 
are similar to the z = 0 sample described in Section |2.2| 
We also characterize the observed and simulated samples by 
star formation, in the form of optical colour. We assign to 
each point in Eigurej^a symbol colour based on the source’s 
U — B colour, where the scaling is identical in the theory and 
data panels. We split the sample at U — B = 1.0 and draw 
separate contours for the blue and red samples. We used 
slightly mismatched Liters in Eigure (rest-frame B or g 
for HST data but rest-frame r for Illustris), because we cre¬ 
ated and analyzed here only the simulated r-band images at 
2 ; = 0.33. In these galaxies, the morphological k-correction 
is small: on average, switching to the B band will increase 
M 20 by ~ 0.1 — 0.2 at most (assuming dust-free), and will 
change G by a negligible amount (< 0.01), compared to the 



_L. 

-10 1 
F (G,M2o) 


Figure 3. Statistical comparison of morphology distributions 
between Illustris (black points and curves) at 2 ; = 0.33 and HST 
data from EGS at 0.2 < 2 ; < 0.4 (blue points and curves). Panel 
a: histogram of F(G,M 2 o). Panels b and c: Position in G-M 20 
space versus F; where the locii differ in Figure these curves 
diverge. Panel d: median absolute deviation of M 20 values at a 
given F, showing how the widths of the main locii vary. Panel 
e: color-morphology plot, showing all original points as well as 
median trendlines versus F. 


r band images. With dust, these differences would be even 
smaller. 

In Figure we examine summary statistics of the EGS 
and Illustris samples. In panel a, the histogram of F{G, M 20 ) 
shows that the two samples both increase in density as F in¬ 
creases above —1, have the same average density to within 
~ 10% at —0.5 < F < 0.5, and then decrease in density 
to nearly zero at F ~ 1. Real galaxies have a more bimodal 
morphology distribution than Illustris galaxies, since the ob¬ 
served number density drops to ~ 2/3 of the peak value at 
F ~ 0.1 where the Illustris distribution is flat. Panels b and 
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Figure 4. Gini-M 2 o as a bulge diagnostic. Left: data from the Zurich Estimator of Structural Types (ZEST; [Scarlata et al.|2007| from 
HST ACS /-band images from the COSMOS survey ( [Scoville et al.|2007| . These data span roughly 0 < 2 ; < 1 and reflect rest-frame g, 
r, and Tband light. ZEST deflned early type galaxies using a principal component decomposition of observed structures, which roughly 
matches the results of visual classiflcation. Right: ~ 40000 mock data points from the Illustris Simulation in the in rest-frame ^f-band, 
for all galaxies with M* > 10 ^-^Mq. Solid red and black lines encircle 68% of the data points, and solid gray lines 95%. /3 is an estimate 
of the bulge fraction derived from the orbits of star particles (Section |4.1| >. This shows how the overall population of Illustris galaxies 
approximately matches the observed morphology distribution. However, the main locus is somewhat wider in Illustris than in the observed 
sample: there are more simulated galaxies with excess G and M 20 than observed. 


c show the position of the G-M 20 morphology locus as a 
function of F : simulated galaxy images have slightly greater 
G and M 20 values than observed ones, as hinted by the Il¬ 
lustris C-M 20 locus in Figure being shifted up and to the 
left. Panel d plots the width of the morphology distribution 
measured by the median absolute deviation from the me¬ 
dian M 20 value at a given F. Outside a narrow range at 
0 < T" < 0.2, the widths of the distributions are almost the 
same with MAD(M 2 o|f) ~ 0.10—0.15. The MAD of the EGS 
sample is large near T" ~ 0.1 owing to the lack of a clearly 
defined locus of galaxies through this region in Figure as 
we noted above, the EGS sample has a much lower number 
density here than in Illustris. Panel e shows U — B color 
versus F for all galaxies in our comparison, with median 
trendlines overplotted. Illustris galaxies are redder by ~ 0.1 
mag at a given F 0 than EGS galaxies, but otherwise the 
extrema are well matched with U — B(F ~ —1) 0.7 and 

U -B{F ^ 1) 1.2. 


In Figure]^ we compare all of our imaged Illustris galax¬ 
ies at 2 ; = 0 with M* > lO^'^M© to the GOSMOS survey 
( Scoville et al.|[20Q7 ) at 0 < z < 1. Here we show results 
for ah galaxies, including both centrals and satellites. The 
~ 10000 model galaxies translates into ~ 40,000 points 
in this plot, with four viewing directions per galaxy. This 
comparison is not perfectly fair, because we are neglecting 
the wavelength dependence of morphology, evolution of the 
galaxy population from z = 1 to z = 0, and sample se¬ 
lection effects. However, since the overall locus in G-M 20 
is similar to the fairer volume-limited sample in Figure 
we conclude that this limitation does not strongly affect the 
general agreement. We show data from Scarlata et al. ( 2QQ7| ), 
who computed morphologies with the Zurich Estimator of 


Structural Types (ZEST), and show the full samples in the 
black contours which represent 68% and 95% of the pop¬ 
ulations. Red contours represent 68% of the “early-types”, 
defined statistically from automated morphologies by ZEST 
and by orbital distribution of stellar mass in Illustris. Thus 
both the real universe and Illustris have similar locations and 
number densities of their late-type and early-type galaxies. 

The main conclusions from Figures and are: 

(i) The overall locus of simulated galaxy points has nearly 
the same size, shape, and location in G-M 20 space as the 
observed one (Figure]^. 

(ii) There exist quenched, early-type galaxies in Illustris 
with a number density relative to star-forming late-type 
galaxies which is within a factor of two from the observed 
one (Figure]^ panel a). 

(iii) Illustris achieves a qualitative match to the color- 
morphology relation (Figures [^, with the extremes of 
the morphology population (\F\ ~ 1) having distinct U — B 
colors within 0.05 magnitudes of observed galaxies, with a 
smooth transition between. 


This is shown, first, by the existence of a large num¬ 
ber of simulated points at M 20 < —2 and G > 0.55, in 
the region occupied by nearby true ellipticals, with very red 
colours (and low SFRs, see Figure and Section 3.2). Sec¬ 
ond, we show contours relative to the same number density, 
and so the fact that the outer blue contours have roughly the 
same position in both panels implies that the relative num¬ 
ber density of bulge-dominated and disc-dominated galaxies 
is reasonably well matched. Finally, Figure [^confirms that 
the number densities as a function of morphology and color 
are broadly matched. 
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These successes are essential first steps that any realistic 
simulation of galaxy formation should satisfy, but which has 
been difficult to attain until recently owing to limitations in 
modeling techniques and in computational resources. This 
agreement implies that Illustris is both large enough to sim¬ 
ulate a fair sample of galaxies, and also realistic enough 
that the number density of early and late types are roughly 
matched to observations. 

We find several basic differences. The theory contours 
in Figure and median statistics of Figure both show a 
tendency to have slightly higher G (by <0.5) at a given F, 
and greater M 20 at a given F (by ~ 0.1), when compared 
with data. As shown by Figure Illustris lacks a clear 
drop in number density (“valley”) between the extremes of 
the morphology and color distributions, causing disagree¬ 
ments of ~ 10 % in number density and ~ 0.1 magnitudes 
in color at a given F ~ 0. This could be a result of our de¬ 
cision to neglect dust (see, e.g., Snyder et al.|2015 ). For the 
simulated galaxies, there are relatively more points in the 
blue contours (late types), reflected by the wider inner blue 
contours. This can also be seen in the simulated mass de¬ 
pendence of quenching ( [Vogelsberger et al.||2014a ). Illustris 
has more “red spirals” of intermediate shape, as reflected 
by the extension of the red contours to the lower left of the 
E/SO/Sa region. 

Despite these differences. Figure [^demonstrates that a 
general correlation between optical morphology and colour 
is in place. Red galaxies tend to be compact spheroids and 
blue galaxies tend to be extended discs. 

3.2 Dependence on star formation 

Observationally, a galaxy’s mass and morphology correlate 
with its star formation rate. Nearby galaxies occur in a bi- 
modal distribution of mass and star formation (or luminos¬ 
ity and colour), in that the luminous/red galaxies appear 
to be a statistically distinct population from star-forming 
galaxies (e.g., Strateva et al.|| 20 ^ ). Surveys have found that 
quiescence is strongly correlated with measures of optical 
structure, such as Sersic index and compactness, albeit with 
large scatter (e.g., Wuyts et al.|2Qlf Bell et al. PM 21 
et al.|2013| |. 


Noeske et al. 

|2007 

Whitaker et al.|2012 ) in the blue squares. 

Torrey et al.[ 

(2014 

) showed that the Illustris model approxi- 


Barro 


Section 13.11 hinted at the overall reasonableness of the 
simulated galaxy shapes and that early-type galaxies are 
redder than late-type ones. In Figurewe plot directly our 
automated non-parametric morphologies from Illustris as a 
function of M*, SFR, and other properties. For the simulated 
galaxies in each bin, we compute the median F(G, M 20 ) us¬ 
ing the definition in Section |2.3| and Figure measuring 
from the unattenuated rest-frame ^-band images. To each 
of these median values of F(G, M 20 ), we assign a colour 
from blue (disc-dominated) to red (bulge-dominated), as 
shown by the legend in Figure We also plot contours 
encircling bins that contain more than 10 measurements, 
black (gray) lines indicating central (all) galaxies. We mea¬ 
sure M* directly from the simulation outputs as the mass 
of the stellar particles contained within a subhalo-centered 
sphere with radius twice the 3D stellar half-mass radius. We 
have checked that this definition does not have a noticeable 
impact on any of our results using M* (especially Section]^. 

In panel (a) of Figure we plot SFR versus M* and 
recover the “main sequence” of star-forming galaxies (e.g.. 


mately reproduces this trend at z = 0 , and we find again that 
its location and slope in this plot is the same as that found 


by Whitaker et al. ( 2012 ). Genel et al. (2014) and Sparre 


et al. (2015) presented the Illustris galaxy main sequence as 


a function of time and show that it roughly matches obser¬ 
vations at 2 : = 0 and 2 ; = 4, but its normalization is too high 
compared with observations at 2 ; ~ 2 . Away from this main 
sequence in Figure at fixed M*, we find that lower-SFR 
galaxies have, on average, earlier structural types as mea¬ 
sured by G-M 2 Q. There are hints that compact or spheroidal 
starbursts populate the highly star-forming sides of this dis¬ 
tribution (see also Sparre et al. poT^ . 

We find a clear separation between bulge-dominated 
and disc-dominated galaxies. While here we are using a 
different morphology diagnostic, the transition from disk- 
dominated to bulge-dominated occurs at a location very 
similar to the one found by Wuyts et al. (2011): at SFR 
~ 0.1 MQyr“^ for ~ lO^^M© and SFR ~ 1 M 0 yr“^ for 
M* ~ lO^^M©. Compared to observations, Illustris appears 
to have too few low-SFR, bulge-dominated central galaxies 
at M* ~ 10 ^°'^M© in panel (a). 

Panel (b) of Figure [^ presents morphology as a func¬ 
tion of overdensity and M*. We use the same definition of 
3-dimensional overdensity 6 as Vogelsberger et al. ( |2014b ), 
who showed that star formation is correlated inversely with 
density at fixed M*, matching observed trends (e.g., Peng 


et al. I polol l. Since F(C, M 20 ) is correlated tightly with 


SFR (panel a), we recover a similar dependence of galaxy 
morphology on mass and environment. Thus, Illustris re¬ 
produces basic features of “mass quenching”, “environment 
quenching”, and the morphology-density relation. 

In panel (c) of Figure we plot SFR/M= 

, a measure of compactness (e.g.. Si .5 from 


versus 


Barro 


et al.|2013 ). This quantity is closely related to a surface mass 


density, which has been shown to correlate with SFR/M= 
out to high redshifts (e.g., Franx et al.|| 200 ^ Omand et al. 


2014). As also seen in individual galaxy tracks by Genel et al. 


(2014), we see that compactness correlates with SFR/M* m 
Illustris central galaxies: the black solid contour in panel (c) 
encircles bins containing more than 10 central galaxies, indi¬ 
cating a tight relationship between and SFR/M*. 

Here, we find that both quantities also correlate tightly with 
bulge strength, as indicated by the change from blue to yel¬ 
low to red bins as one moves along the relation to the upper 
right corner. Gompared with results by Barro et al. I poT^ at 
2 : > 0.5, Illustris galaxies have somewhat smaller SFR/M* 
and This difference owes to the decline of SFR/M* 

over time (e.g., Noeske et al.|[2007 [Genel et al.||2014| and 
larger-than-observed Illustris galaxies (Section [5.1| . 

Another feature of panel (c) in Figure [^ is the struc¬ 
tural difference between central and satellite galaxies, which 
can be seen by the difference between the gray and black 
contours. These encircle all and central galaxies, respec¬ 
tively. There is a single tight correlation between SFR/M* 
and M*F“^ 2 ^ in central galaxies, but less so among satellite 
galaxies, which are quenched and early type but have lower 
at a given SFR/M*. 

With the three trends of Figure [^ we conclude that 
the Illustris simulation produces a 2 : = 0 galaxy population 
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Figure 5. SFR versus mass (panel a, top left), overdensity versus mass (panel b, bottom left), and SFR/M* versus compactness (panel 
c, top right), with colours proportional to galaxy structure (Equation]^. In each bin, we measure the median F(G, M 20 ) value from 
rest-frame g-band images and assign it a colour following panel d. Bulge-dominated galaxies are colour-coded red and disc-dominated 
ones are blue. The gray (black) contour outlines a region in which each bin contains 10 or more (central) galaxies. The green solid line 
in panel (a) is the star-forming main sequence parameterized by [Whitaker et Compared to surveys, Illustris recovers roughly 

the s ame dependence of structural type o n and SFR (e . g., [K auflhna^ n et al. [|20Q3[ [Wuyts et al.||2011| >, on 5 and M* ( |Peng et al.| 

Franx et al. 2008 Barro et al. 12013 Omand et al.l2014 1. These trends appear to be a natural 


20101, and on SFR/M* and (e.g 

consequence of galaxy formation processes 


simulated with physics models crafted to match global star formation histories and stellar 


mass functions. 


with optical morphologies that successfully match trends 
emphasized in survey analyses. These include the SFR-M*- 
Morphology relation (e.g., Kauffmann et al^|2003 Blanton 
et al. 2Q03| Wuyts et al. 2011), environment versus mass 
quenching (e.g., Peng et al.|2010 ), and the tight correlation 
between galaxy compactness, star formation, and morphol¬ 
ogy ( e.g., [Franx et al.|2QQ8 Williams et al.|20lQ Barro et al. 


2013 Oniand et al.||2014 ). These trends have been thought 

important for understanding key issues in the formation of 
massive galaxies, such as bulge formation, quenching, and 
emergence of the Hubble Sequence. 


The recovery of the median morphology of simulated 
galaxies as functions of these properties, trends not explic¬ 
itly matched in the Illustris physics model selection, is an 
intriguing success of such simulations. The ingredients re¬ 


quired to obtain reasonable global stellar mass functions and 
star formation histories appear to also imply these more de¬ 
tailed relationships between galaxy morphology, mass, and 
star formation. 


3.3 Dependence on stellar mass 

In this section we show in more detail how the morpholo¬ 
gies of Illustris simulation galaxies depend on stellar mass. 
Figure [^presents the distribution of G-M 20 in Illustris as a 
function of stellar mass at z = 0, as measured from the rest- 
frame g band images. We separate the measurements into 
three bins from 10.0 < log^Q M^IMq < 12.0. In each panel, 
we bin galaxies by G-M 20 values and show an example im¬ 
age selected from each bin as having most nearly the median 
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Figure 6. Rest-frame ^f-band morphology distributions and colour-composite images as a function of mass. For each bin, we show the 
galaxy having F(G, M 20 ) nearest the median value. White (gray) contours encircle bins containing 100 (10) or more galaxies. From these 
we see an expected trend: at lower mass (bottom), most galaxies are still forming stars in a disk, and therefore reside in the late-type 
region (M 20 > —2, G < 0.6). At the highest masses (top), galaxies are quenched, and therefore they are concentrated mainly in the 
early-type region. However, at log^^Q M*/Mq < 11.0, we identify a peculiar population of galaxies: there are two peaks in G-M 20 in the 
middle panel, one at the expected location bridging the early/late type divide as galaxies evolve and quench, and another off to the left, 
at M 20 > —1.0 and G ~ 0.4. See the text for more discussion. 
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morphology, F{G^ M 20 ). White (gray) contours encircle bins 
containing more than 100 (10) galaxies. 

The basic trend of Figure [^follows what we would ex¬ 
pect given that morphology in Illustris, and most galaxy 
formation models, is closely correlated with mass. Simulated 
galaxies tend to occupy the Sb/Sc/Irr region of the G-M 20 
diagram at low mass, and the main locus shifts smoothly to 
the E/SO/Sa region at M* > 10^^ M©. This overall shift is 
visible in the galaxy images: at the top are quenched bulge- 
dominated galaxies, at the bottom star-forming galaxies. 
The star-forming galaxies have a diverse morphology: some 
are almost entirely blue and disc-dominated, while others 
are composites, having a compact red bulge. 

At first glance, there are two possible issues with the 
Illustris model in G-M 20 versus mass: 


(i) It appears there are very few, if any, galaxies with 
M* < Mq with a bulge-dominated morphology 

(M 20 < —2, G > 0.55); compare also our Figure panel 
a to Wuyts et al.| ( |2011| . We believe that this relates to 


the difficulty of properly quenching star formation in some 
galaxies with M* < 10^^ M©. 

(ii) At < M*/Mo < 10^^, the G-M 20 locus appears 

to separate into two distinct populations. The main one near 
the center of the diagram, intermediate between the adjacent 
panels, and another with very extended morphology: M 20 > 
-1, G < 0.45. 


The origin of the second locus of galaxies with masses 
M* < 10^is unclear. The feature appears to be caused 
by rings of star formation which are visible in the example 
images in Figure]^ (also Fig ure p^ . We discuss the origin of 
this phenomenon in Section 

It is easy to see how M 20 can be sensitive to this feature. 
In short, M 20 represents the spatial extent of the brightest 
quintile of a galaxy’s pixels. If there exists a ring of star 
formation that is brighter than in typical spirals, as we hy¬ 
pothesize here, then the brightest (star-forming) pixels may 
be unusually separated on the sky. Figure shows an abun¬ 
dance of objects of this character in the rest-frame g-band, 
but this feature is also present in the rest-frame i and H- 
band morphologies and images. The images highlighted from 
this region show both blue pixels and yellow/orange pixels, 
indicating that these “rings” can be sustained long enough 
to create older stellar populations at that radius. 


4 OPTICAL STRUCTURE AND KINEMATICS 


The optical morphology of an image is a projection of the 
SEDs of stellar populations orbiting in the galaxy. Since sim¬ 
ulations predict these populations, we can relate structural 
parameters of simulated galaxy images to the motions of 
stars. In principle, the kinematic information provided by 
large-scale simulations could be used not only to interpret 
observations (e.g., Kassin et al.||2014 ) but also to constrain 
galaxy physics models from integral held unit surveys. More¬ 
over, optical and kinematic morphologies are sensitive in dif¬ 
ferent ways to the stellar mass and SFR distributions, and 
understanding this connection is important for large imag¬ 
ing surveys. 


4.1 Morphology versus Rotation 


We hrst compute several measures of galaxy rotation based 
on stellar motions. We follow previous studies (e.g., Scanna- 


pieco et al.|2010 Sales et al.|2012 Marinacci et al.|2013 ) to 

parameterize the orbits of star particles, and contrast these 
quantities with photometric morphologies. 

For each galaxy subhalo, we compute the total specihc 
angular momentum j^ — ^ (rn v x r)/^moi the star par¬ 
ticles within ten times the galaxy’s stellar half-mass radius. 
This radius is chosen to correspond with observational stud¬ 
ies (e.g.. Fall Romanowsky||2013 ), which hnd that a sig- 
nihcant fraction of a galaxy’s stellar angular momentum lies 
beyond the half-mass radius ( Romanowsky Fall|[20T2 ). 
Genel et al. (2015) found that using all star particles in 


the subhalo instead of restricting to IOR 1/2 (or 5 R 1 / 2 ) does 
not change the resulting scaling relations in Illustris galax¬ 
ies. We define the direction of to be the 2 : unit vector 
and compute the specific angular momentum of each star 
particle in that direction, jz = (v x r)^. We then define a 
circularity parameter e = jzljz{E)^ where E is the particle’s 
binding energy, and jz{E) is the maximum specific angular 
momentum among the 100 particles with a binding energy 
closest to E. Finally, for each galaxy subhalo we compute 

= |j*l several summary statistics of the distribution 
of e. These quantities correlate with rotational support or 
bulge mass fraction. 

In Figure we plot F(G, M 20 ) of massive galaxies 
(M* > 10^^M©) against three such parameters: 

(i) [3 — twice the fraction of star particles with e < 0 (e.g., 
Abadi et al.||2003 ). Perfectly cold discs should have ~ 0, 
and bulge-dominated systems have /3 ~ 1. 

(ii) login. 7 * Since specific angular momentum 

scales with (Romanowsky & Fall 2012), j^/MV^ is 

a measure of relative rotational support. This quantity is 
similar to projected, observed par ameters such as \r used 
in kinematic surveys (e.g., Atlas^^ Cappellari et al. 2011b), 
and has units of km s“^ kpc 

(hi) /^rot = defined by Sales et al. 

(2012), where K is the total kinetic energy of star particles 


and R is the distance between the star particle and galaxy 
center. Krot is the fraction of kinetic energy occupied by 
ordered rotational motion. Bulge-dominated systems have 
^rot ^ 1, while perfectly cold discs have hirot ~ 1- 


In each case, E(G,M 2 o) has some dependence on the 
kinematic measure. Simulated galaxies with /3 > 0.5 are 
optically early type {E{G, M 20 ) > 0.0), and there is a 
locus of objects with ~ 0.2 dex scatter extending from 
{/3 - 0.3, F{G,M 2 o) - 0) to {13 - 1, F{G,M 2 o) - 0.7), 
so that F(G, M 20 ) correlates positively with j3. There is a 
large diversity of optical morphologies when /3 < 0.5, with 
F{G,M 2 o) spanning from —1 to 0.5. 

Similarly, the rotation measures and Krot 

are inversely proportional to optical diagnostics of bulge 
strength. F(G, M 20 ) versus Krot roughly follows the quan¬ 
tity I — j3, with a wide scatter in disc-dominated galaxies, 
albeit slightly less scatter than for /3. This scatter appears 
to increase in galaxies with lower mass (M<10^°'^M©), and 
could reflect more complex structures in these galaxies com¬ 
pared with very massive ones. For example, the [3 indica¬ 
tor assumes that bulges are rot at ion-free, but we know that 
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Figure 7. Simulated distribution of optical morphology versus kinematic bulge fraction (left), ratio of specific angular momentum to 
(center), and ^Crot parameter (right), in galaxies with M* > 10^^ Mq. White (gray) contours enclose regions containing 70% 
(95%) of the model sources, and colours of the rectangles indicate the fraction of quiesc ent galaxies in each bin. We define quiescent 
fraction as the fraction of galaxies with log^^Q SSFR/(yr“^) < —0.24(log]^Q M*/Mq) — 8.50 (j Oinand et al.|2014| . We define (3 as twice the 


fraction of star particles with negative circularity parameter; see text for full definition. Section 4.1 defines and Urot- Light-based 

morphology traces kinematic-based ones for > 10 ^(right), especially at /3 > 0.4. Since the distribution of colours is largely 
horizontal, the optical morphology traced by G-M 20 is tightly correlated with quiescence. 


many early-type galaxies are significant rotators. Addition¬ 
ally, structures like bars and counter-rotating disks can con¬ 
found the kinematic indicators used here. 


4.2 Quiescence versus Morphology and Rotation 

In Figure we colour bins by quiescent fraction, de¬ 
fined as the fraction of galaxies with log^Q SSFR/(yr“^) 
< —0.24(logiQ M^^/M©) — 8.50 (e.g., Omand et al.|[^14| ). 
This roughly matches the definition of quiescent or passive 
galaxies used in several other studies (e.g., Brammer et al. 
Woo et al. 2015| . As we showed in Figure opti- 


2009 


cal morphology traced by F(G, M 20 ) correlates tightly with 
quiescence at almost any mass M* > 10^°M©. 

Intriguingly, quenching appears to be correlated less 
well with kinematic morphology: the bin colours are dis¬ 
tributed horizontally in the T(G, M 20 ) versus /3 and /^rot 
panels of Figure By contrast, the colours trace a more 
diagonal trend in F versus this x-axis quantity 

helps to predict quenching at fixed F. However, this is true 
primarily in rotation-dominated systems. 

This implies that optical morphology closely traces star 
formation activity in Illustris, even when the simulated kine¬ 
matic morphology does not. This could result from several 
factors, including that the kinematic tracers are mass-based 
while the image tracers are light-based. Therefore the op¬ 
tical profile follows the mass profile only if star formation 
does. 

Figure [^summarizes how F(G^ M 20 ), SFR/M*, and ro¬ 
tation correlate with each other in Illustris. In the left panel, 
we see that at fixed relative rotation, F is inversely propor¬ 
tional to SFR/M*. And at fixed SFR/M*, F is inversely 
proportional to In the right panel, we over-plot 

images chosen to have most nearly the median F value of 
galaxies in the associated bin. 


In Figure the low-redshift Illustris galaxies occupy a 
sequence very similar to the one put forth by observational 
studies (e.g., [Emsellem et al. 2007| Cappellari et al. 2007 


2011b): quenched bulge-dominated galaxies with very little 


rotation are rare and occupy the bottom left tip of the locus 
of galaxies. Moving right from these slow rotators we see 
the “fast rotators”, which can be visually very similar to 
the slow rotators when viewed face-on. However, there are 
obviously some disc-dominated galaxies viewed edge-on in 
the lower right quadrant of Figure as in real galaxies. 
In rotating galaxies, there is a continuum of star formation 
activity which is correlated with the appearance of spiral 
arms and with having a more disc-dominated distribution 
of optical light. 

The trends in this figure follow those which kinematic 
surveys of nearby galaxies have identified recently. For ex¬ 
ample, galaxies with bulge-dominated F values occupy a 
huge range of relative levels of rotation. Thus, in red galax¬ 
ies, it is very difficult to tell from optical morphology alone, 
either visually or automatically, whether a galaxy has sig¬ 
nificant rotation. Moreover, at fixed nonzero rotation level, 
simulated galaxies have a wide range of SFR activity. At 
minimal rotation (log^Q + 5 < 0), most simulated 

galaxies are quenched. This is reminiscent of recent work em¬ 
phasizing a kinematic perspective on the Hubble sequence 


(e.g., Emsellem et al.|20Q7 Krajnovic et al.|2008| Cappellari 


|et al. 2011b), whereby early type galaxies with high levels of 

rotation are more closely related with “anemic spirals” (e.g., 
van den Bergh|[T976 ) than to true slowly-rotating elliptical 


galaxies. 


4.3 Morphology-density 


Vogelsberger et al. (2014a) showed how quiescence traces 


environment in Illustris galaxies, and we showed in Figure 
that this translates to the expected dependence of galaxy 
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Figure 8. SFR/M* versus rotation in Illustris. Left: We colour each bin by the median F(G, M 20 ) of the galaxies contained therein. 
Right: We fill each bin with the colour-composite image having most nearly the bin’s median F(G, M 20 ) value. Therefore, we obtain a 
kinematically defined galaxy sequence qualitatively matching observed ones (e.g., [Cappellari et al.||2Qllb| >. At the bottom left are the 
slowly rotating early-type galaxies. Moving to the right, the automated or visual classifications remain early-type, while the rotational 
measure increases and several edge-on examples appear disc-dominated. Having high rotation and high SFR/M*, the top right of the 
figure contains obviously late-type star-forming galaxies as well as the ring galaxies of Section |3.3| 


structure on density. In Figure we explore in greater de¬ 
tail such a “morphology-density” relation, and how it ap¬ 
pears in Illustris when using the slightly finer definitions of 
morphology discussed in Section |4.2[ We use p as the vol¬ 
ume number density of galaxies following the conventions of 
observational studies (e.g., Dressler||1980| [Cappellari et al. 


2011b). 


First, in the top panel of Figure we classify simu¬ 
lated galaxies as late-type discs, early-type discs, or ellipti¬ 
cals using galaxies’ SFR/M* and F(G, M 20 ). At low densi¬ 
ties, p < 1 Mpc“^, the population is > 80% star-forming 
disc-like galaxies {F < 0.3, SFR/M* > 0.01 Gyr“^). As 
the density increases, the share of late-type discs shrinks 
smoothly, falling below 0.5 when p > 10 Mpc“^, while the 
share of early-type discs smoothly increases. The share of 
near or pure elliptical galaxies {F > 0.3) is small but in¬ 
creases smoothly to ~ 10% over this range. 

To choose these classification boundaries, we followed 
an approach similar to those taken by kinematic studies, 
which have highlighted the difference between slowly rotat¬ 
ing ellipticals and red, rotating disc galaxies (e.g., Cappellari 


et al. 2011b), both of which could be classified visually or 


automatically in projected light as early type. Galaxies with 
F > 0.3 are almost entirely classical bulge-dominated ellip¬ 
ticals, while those with F < 0.3 are a mix of bulge- and disc- 
dominated galaxies, often with significant rotation. 

In the bottom panel of Figure we classify galaxies 
directly using their stellar mass kinematics and SFR/M*. 
We follow an approach similar to above, selecting slow rota¬ 
tors as those with extremely low specific angular momentum 
relative to their mass: log^Q + 5 < 0.0. These are 

almost exclusively quenched, massive, and optically early 
type {F > 0.0: Figure]^. By contrast, as we saw in Sec¬ 
tion |4.1| simulated rotating galaxies have a broad range 
of properties. Therefore, we separate the rotating galaxies 


at SFR/M=, = 0.01 Gyr“^ into quenched and star-forming 
classes. 

We find that the kinematic morphology-density relation 
is very similar to the optically defined one: At p < 1 Mpc“^, 
the population is > 80% rotating star-forming galaxies. As 
the density increases, the share of rotating quenched galax¬ 
ies smoothly increases, reaching 50% at p ~ 10 Mpc“^. The 
share of slowly rotating galaxies increases smoothly from 
~ 1% to ~ 10% over this range. At the lowest densities, 
there is a hint that the simulated galaxy population is satu¬ 
rating as star-forming discs, in contrast to observations (e.g., 
Gappehari et al.||2011b ). However, this is an expected con¬ 


sequence of the fact that the Illustris model yields too little 
quenching in the low-mass central galaxies dominating these 
fractions. 

In both panels of Figure the share of early-type 
discs equals the share of late-type discs at a volume den¬ 
sity p ~ 10 — 20 Mpc~^, closely matching the results of 
Atlas^^( Gappehari et aL]|2011b ), which has a very similar 
selection limit (M* > 6 x 10^ M©) to the Illustris sample here 
(M* > 4.5 X 10^ M©). The transition from star-forming (low 
densities) to quenched (high densities) is extremely smooth, 
with no evidence for a break or dramatic change at any 
density threshold, even beyond the Atlas^^ sample range 
at logj^Q p> 1.5. Moreover, the fraction of slow rotators, rel¬ 
ative to ah quenched galaxies, is nearly constant at a lmost 
20%, also roughly matching the results of Atlas^^ (Gap- 
pehari et al.||20lTb ). The normalizations of both simulated 


and observed ratios is somewhat arbitrary, since they de¬ 
pend on the selection cut between slow and fast rotation. 
Regardless, we find that the simulated iVsiow/Wast is nearly 
constant over three orders of magnitude in density. 

Above a density of p ~ 100 Mpc“^, the statistical un¬ 
certainty on the population fractions is too large to continue 
following these trends, owing to the finite volume of Illustris 
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Figure 9. Simulated morphology-density relations. Top: Type 
fractions, classified by SFR/M^c and F(G, M 20 ), as a function 
of local galaxy volume number density p for M* > 10 ^-^Mq. 
Bottom: Type fractions, classified using SFR/M* and rot ation. 


We follow integral field unit surveys, such as Atlas^^ (Cap- 
[pellari et al.|[^lla| >, to classify simulated galaxies as slow rota¬ 
tors, rotators with star formation (a.k.a. late-type galaxies), and 
quenched rotators. The share of late-type discs falls smoothly as 
density increases, falling below the share of early-type discs at 
p 10 — 20 Mpc“^. We plot the ratio between the numbers 
of slow and fast rotators among quenched galaxies (slow rota¬ 
tors are almost exclusively quenched, i.e. Figure]^, from Illustris 
(black line with error bars indicating the Icr Poisson noise) and 
the global average from |Cappellari et al.| ( [2011b| dotted line). At 
much higher densities (i.e., the richest clusters), there are insuffi¬ 
cient statistics in the z = 0 Illustris snapshot to continue following 
this relation. 


(106 Mpc)^. While this is adequate for comparing to very 
local volumes (e.g., Virgo cluster samples), we require sim¬ 
ulations with volumes > (200 Mpc)^ to probe this regime 
more with greater statistical significance and contrast with 
data from forthcoming wide-area, high-resolution imaging 


surveys from (e.g., Euclid and WFIRST, Laureijs et al. 
WI]|Spergel et al.|2015| |. 


However, to improve the situation at lower galaxy 
masses (and lower number densities) may require qualita¬ 
tively improved quenching mechanisms. While the colours of 


Illustris satellite galaxies are in agreement with data (Sales 
eTaTIpOMl l, it is possible that numerical resolution better 


than ~ 1 kpc is required to implement efficient enough feed¬ 
back physics, or to allow the ISM to quench through dif¬ 
ferent mechanisms in isolated or central galaxies (See also 
Section |6.4|) . 


5 MORPHOLOGY, QUENCHING, AND 

GALAXY PROPERTIES 

In Sections and ^ we showed how optical galaxy mor¬ 
phologies in the Illustris simulation at z = 0 are a reasonable 
match to several observed trends broadly related to massive 
galaxy evolution. This appears to be a consequence of craft¬ 
ing the sub-grid supernova and AGN feedback models to 
roughly recover the total stellar mass content of halos over 
cosmic time. An interesting question becomes: Given such 
tuning, is the success with galaxy morphology inevitable, or 
is this particular class of comprehensive models special in 
some way? 

In this section, we explore how Illustris galaxy mor¬ 
phologies depend on other aspects of the simulated galaxies, 
such as their optical sizes (Section |5.I| ), SMBH masses (Sec¬ 
tion 5.2), and dark matter (DM) halo masses (Section |5.3|). 


The idea is to search for predictions that could be unique to 
the specific galaxy physics assumed in Illustris. 

Figure presents several correlations among galaxy 
properties. In addition to raw correlations, we over-plot con¬ 
tours indicating the “quiescent fraction” of galaxies in each 
of the 2-dimensional bins, in order to give some idea of the 
relationship between these quantities and galaxy star forma¬ 
tion (i.e.. Section 3.2). We define quiescent fraction as the 
fraction of galaxies satisfying: 

logio SSFR/(yr-') < -0.24(logio M./M©) - 8.50, (7) 

as used also in Section 4.2 and following Omand et al. ( 20I4| ). 
This corresponds to SSFR ~ I0“^^yr“^ in galaxies with 
M* ~ lO^^M©. All sizes and morphologies presented in this 
section are measured in the rest-frame SDSS-^ band images. 


5.1 Size-mass-morphology 


We show how size, mass, and quenching are related in the 
bottom left panel of Figure Here we make a rough com¬ 
parison to the distribution of SDSS galaxies presented by 
Omand et al. ( 20I4| ). Genel et al. (in prep.) will study the 
size distributions in more detail. We use sizes measured from 
the rest-frame g band images, as described in Section |2.3| 
There are thus four sources with morphology and size mea¬ 
surements for each simulated galaxy, one each for the four 
default random viewing angles. 

The Illustris size-mass correlation is shallower than the 
observed one for M^/M© > 10^'^. At M*/M© ~ 10^^, the 
observed and simulated distributions are very similar. How¬ 
ever, observed galaxies shrink more quickly as M* decreases 
than do Illustris simulation galaxies. The lower end of the 
size distribution at R 1/2 ~ I kpc may be affected by the rel¬ 
atively coarse spatial resolution (^0.7 kpc) of the simulated 
galaxies. Indeed it may be impossible for Illustris galaxies to 
reflect a size smaller than a few times this resolution scale. 
However, this effect alone would not necessarily make the 
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Figure 10. Simulated relationships between size, mass, morphol¬ 
ogy, and star formation. Sizes and morphologies are as measured 
from the unattenuated rest-frame g band images. White (gray) 
contours encircle bins with 50 (10) or more simulated points, and 
the rectangle colours reflect the fraction of quiescent galaxies in 
each bin. In the lower left panel, we plot half-light radius versus 
stellar mass, showing that the simulated distribution is similar to 
the observed one at M* ~ 10 ^ ^Mq but too large by a fa ctor ^ 2 
at IO^^Mq. Illustris recovers the observed trend (e.g., |Omand| 
let al.|[M4t that the distribution of quiescence is tilted with re¬ 
spect to the M* axis. This means at fixed M*, smaller galaxies 
are more likely to be quenched, known to be true for both nearby 
( |Kanffmann et al.|2Q03| > and distant ( |Barro et al.|2Ql^ galaxies. 
The middle row reflects that morphology is also an effective pre¬ 
dictor of quenching, both at flxed (left) and also at flxed i^i /2 
(right). The top panel plots these same variables against SMBH 
mass. We And that objects residing on or above the mean Mbh" 
M* relation are slightly more likely to be quenched (top left), as 
expected and found by |Sijacki et al.| ( |2Q15| >. 


largest simulated galaxies significantly larger. Thus we con¬ 
clude that the upper end of the simulated galaxy R 1/2 values 
are a factor of ~ 2 larger than observed at M^jM^ ~ 10^°. 

Although Illustris has fewer quenched/compact galax¬ 
ies (or larger galaxies on average) at M* < the 

diagonal shape of the quenching contours in the M^-Ri /2 
plane is nearly a perfect match to the results of |Omand| 
et al.| ( [20T4 ), where quenching is most tightly correlated with 


3.2 


M^lRil 2 with a ~ 1-2. As we showed also in Section 
this implies that quenching tends to trace the compactness 
of stellar density profiles in roughly the correct manner in 
the Illustris simulation. 

The left middle panel of Figur e shows G-M 20 ver¬ 
sus mass, also discussed in Section |3.3| The Illustris sim¬ 
ulated galaxies reflect the well known rough correlation of 
automated morphology with stellar mass or luminosity (e.g., 
Kauffmann et aL||2003 Strateva et al.[|20Ql| ), with morpho¬ 


logically early type “red sequence” at M^/M© > 10^^. As in 
the Ri /2 — AT* panel, horizontal or diagonal contours in the 
G-M 20 —M* panel imply that morphology is an informative 
predictor of star formation at fixed M*. In fact, the cen¬ 
tral panel of Figure implies that a galaxy’s morphology 
is also a better predictor of quenching in Illustris than its 
radius alone. At fixed R 1 / 2 , average star formation activity 
is inversely proportional to F(G,M 2 o). 


5.2 SMBH-morphology 


Since the AGN feedback model suppresses the low-redshift 
SFR of galaxies, we might expect quantities correlated with 
the SMBH output energy (for example, their mass) to be 
also correlated with star formation activity or quenching. 
Indeed, Sijacki et al. ( |2015| demonstrated that at a given 
M*, galaxy colour (SFR) is directly (inversely) proportional 
to Mbh- In the top row of Figure we show the predicted 
relationship between SMBH mass (Mbh) and optical mor¬ 
phology. 

In the top left panel, we recover the same result as|Si-| 


jacki et al. (2015): at fixed M*, galaxies with higher Mbh 


are more likely to be quiescent. There is also a hint that 
at fixed Mbh, lower M* implies more quiescence, a feature 
that we explore in more detail in Section [5.3| Similarly, at 
fixed Ri/ 2 , and at fixed F{G, M 2 o)^ 0, galaxies with higher 
Mbh are more likely to be quiescent (top center and right). 
However, at intermediate morphology {F ~ 0), and at any 
Ri/ 2 , at fixed Mbh, higher F or lower R 1/2 implies more 
quiescence, on average. 


5.3 Halo mass-morphology 


From single trends alone, it is difficult to assess directly the 
impact of astrophysical processes on the resulting galaxy 
structures. In part, this is because these processes alter the 
structural quantities themselves, and so the baryons in a 
given galaxy may have a very different history under differ¬ 
ent model assumptions. However, a galaxy’s total mass is 
likely far less sensitive to this issue. 

Therefore, investigating how closely quenching and 
morphology depend on halo mass in Illustris may be a pro¬ 
ductive avenue for communicating its testable predictions. 
For example, if the mechanism most responsible for setting 
galaxy morphology depends purely on DM halo mass (Mh), 
then we might expect there to be no residual correlations 
among important galaxy properties at fixed Mh- In con¬ 
trast, some models assume that baryonic processes such as 
feedback regulate galaxy evolution, while their implementa¬ 
tion may also depend on the galaxy’s halo properties (e.g., 
Vogelsberger et al. 2014b). Therefore it is uncertain how 
halo mass might relate to galaxy morphology in a simulation 
like Illustris. Even assuming that all comprehensive models 
matching global stellar mass quantities also reproduce basic 
morphology trends (Section |^, galaxy properties at fixed 
halo mass might be an effective means by which to choose 
among them. 


In Figure 11 we plot E(G, M 20 ) and compactness 
(M*R“^ 2 ^; Section 3.2) against total mass Mh and quies¬ 
cent fraction, as defined in equation 0 - For a given galaxy, 
we define Mh as the total mass of all particles (all types) 
and cells bound to the subhalo, not including mass associ¬ 
ated with subhalos of this subhalo. Qualitatively, Illustris 
galaxy morphologies depend on Mh in a very similar way 
to observed low-redshift galaxies (e.g.. Woo et M^|2015 ): at 
fixed Mh, optical structure is very strongly correlated with 
whether a galaxy is quenched or not. However, it is difficult 
to immediately assign causality. 

Over time, the existence of quenching implies that a 
galaxy will form fewer stars than it would have without the 
quenching mechanism(s). This implies that in a given halo. 
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Figure 11. Quiescent fraction and optical morphology versus 
total mass for massive galaxies in Illustris, with colours rep¬ 
resenting the fraction of quiescent galaxies in each bin, as defined 
in equation 0 At a given M/j,, the quiescent fraction increases 
(bins become red) as F or increases. Illustris galaxies 

at fixed Mfi are much more likely to be quenched if they have 
a dominant central bulge component, as inferred in real galaxies 
(e.g., [Woo et al.|2015| >. 


quantities which trace the integrated star formation history 
(SFH) could provide a direct constraint on the key physical 
processes. For example. Figure shows how morphology 
depends on M* in simulated massive central galaxies at fixed 
Mh . The bin colours as a function of Mh and M* (top panel) 
indicate that at Mh ~ IO^^Mq, galaxies with lower M* have 
higher F{G, M 20 ), on average. In other words, galaxies with 
higher M* at fixed Mh are more disc-dominated. |Pillepich| 
et al. pMIl i also discussed this idea, in which these more 


massive disc-dominated galaxies also have shallower stellar 
halo mass profiles than their less massive, bulge-dominated 
counterparts at fixed Mh- 

The bottom panel of Figure shows explicitly the 
residual correlation between F and M* at fixed Mh- To de¬ 
fine this relation, we fit the F versus log^^Q Mh and log^^Q M* 
versus log^Q Mh correlations with orthogonal linear regres¬ 
sion (Boggs & Rogers 1990) and subtract the predicted 
from actual values: AF = F — F|m^, AM* = log^Q M* — 
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Figure 12. Morphology versus stellar mass at fixed halo mass for 
central galaxies in Illustris. Top: relation, with bin colours 

representing optical morphology, F(G,M 2 o). Bottom: Residual 
correlation between F(G, M 20 ) and M*, holding Mh fixed; see 
text for definitions. In this panel we colour the bins by quiescent 
fraction. Black lines are orthogonal distance regression ( [Boggs ^ 
[Rogers! 1990| > fits to the raw data. We find a statistically signifi¬ 
cant (log^^Qp <C —3) negative residual correlation between F and 
M* at a given Mh , implying that galaxies with higher stellar mass 
have later morphological type, on average. We interpret this as 
a consequence of the Illustris physics model suppressing galaxy 
star formation: in a given halo, less quenching leads to later type 
galaxies and more total stars. More quenching leads to earlier 
types and fewer total stars. 


logj^QM^lM^. With these, we can understand the internal 
correlations among galaxy properties and determine which, 
if any, predict galaxy morphology better than others (e.g., 
Hopkins et al.pOO? ). For the fitting procedure, we assume 


uniform measurement uncertainties of 0.1 dex for masses 
and 0.2 for F, so that the uncertainties are a similar frac¬ 
tion of the intrinsic scatter in each quantity. For clarity, we 
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show AF versus AM* only for galaxies in the halo mass 
range 12.0 < log^QM/^/M© < 12.5, which is a range dur¬ 
ing which the simulated galaxy population transitions from 
mostly star-forming to largely quenched, and therefore we 
might expect the most obvious consequences. 

We find a negative residual correlation between M* 
and F at fixed Mu in 10^^ < MhIMq < 10^^'^. Ac¬ 
cording to t-tests, the null hypothesis (no correlation) is 
overwhelmingly unlikely to yield these data. We also split 
the sample into star-forming and passive subsamples, at 
logio(S^^^/^*) — “2.0, and find in each a negative cor¬ 
relation at high significance. 

Therefore, in halos of the same total mass, simulated 
galaxies with above-average stellar mass are more disc- 
dominated than average, and galaxies with below-average 
stellar mass are more bulge-dominated than average. The 
reason could be simple: stars typically form in a disc config¬ 
uration, and therefore if there is more total star formation in 
a particular galaxy, there will be both more stars and their 
morphology will be more disc-like, on average. 

The colours in the bottom panel of Figure indicate 
quiescent fraction, which tends to be higher in galaxies that 
are more bulge-dominated than average (higher AT). In 
contrast, the average current SFR/M* or quiescent fraction 
is roughly constant as a function of AM*. This was also 
seen in the red fractions of massive simulated galaxies by 


Vogelsberger et al. (2014b) and quenched fractions in SDSS 


by Woo et al. (2012). Therefore, galaxy morphology may be 
a clearer diagnostic of the time-integrated effect of feedback 
in Illustris than current SFR. 

A candidate for the cause of this effect is the radio-mode 
AGN feedback, which is the mechanism most responsible 
for reducing star formation at late times in massive Illustris 
galaxies (Vogelsberger et al. 2013). If so, we might expect 
galaxy properties to have residual correlations with the sim¬ 
ulated SMBH mass at fixed Mh- Indeed, Sijacki et al. (2015) 
found that galaxies with over-massive SMBHs are redder 
than average, which we showed also in Figure with re¬ 
spect to M*. 

In Figure |13| we present the MBH-M/^-morphology re¬ 
lation in Illustris. The top panel shows that above-average 
Mbh implies higher median T, as expected from any rea¬ 
sonable SMBH-bulge correlation. Since F and M* are neg¬ 
atively correlated at fixed Mh (Figure p^, it is likely that 
these same galaxies with above-average Mbh have below- 
average M*. The bottom panel of Figure [^confirms this: at 
fixed M/j, simulated galaxies with above-average Mbh have 
below-average M*. In addition, above-average Mbh implies 
earlier-than-average morphological type (redder bins). 

If true, then the effect of AGN feedback on galaxy mor¬ 
phologies at fixed Mh is indirect. It results from the two 
facts that 1) AGN feedback quenches star formation over 
long periods, and 2) star formation over an extended period 
leads to disc-dominated morphologies. Differences in evolu¬ 
tion of baryons in similar-mass halos would lead to a varying 
integrated feedback efficiency. This would have an effect in 
the sense described above for the M*-M/j-morphology rela¬ 
tion. Since morphology is so tightly linked with the extended 
star formation history (SFH), other physical processes act 
less signficantly on morphology, and the overall SFH poten¬ 
tially masks these effects. For example, a merger may cause 
a late-type galaxy to temporarily evolve into an early-type. 
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Figure 13. Stellar mass versus black hole mass at fixed halo 
mass in Illustris, following same procedure as Figure |12| Top: 
Mbh-M/j relation. Bottom: Residuals of M* and Mbh at fixed 
Mh- Bin colours represent optical morphology, F(G^ M 20 ). We 
find a statistically significant (log^^Q p <C —3) negative correlation 
between AM* and AMbh, implying galaxies with higher SMBH 
mass have lower stellar mass. Again, we can interpret this as a 
consequence of the Illustris physics model: In a given halo, higher 
SMBH mass implies more total energy deposited by feedback and 
a greater reduction in the total stellar mass formed. 



However, if quenching isn’t sufficient owing to an undermas- 
sive SMBH, then this galaxy might continue forming stars 
and re-grow a disc over a few Gyr (c.f., Snyder et al.|2015 ), 
leading to a higher M* and more disc-dominated morphol¬ 
ogy at late times. 


5.4 M* and quiescent fraction 

If we assume that quenching determines M* at fixed Mh 
in Illustris in the sense described in Section [5ffil one conse¬ 
quence is that we should be cautious with quantities depend- 
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Figure 14. Residuals of compactness and M* at fixed halo mass, 
with bin colours representing quiescent fraction. cor¬ 

relates positively with AM*, in direct contrast to the behavior of 
AF in Figure 



ing directly on M*. For example, the compactness quantity 
we used above correlates positively with M* at fixed in 
contrast with F(G, M 20 ) and Mbh; see Figure 14 Thus, it 
is possible that using quantities directly proportional to M* 
can obscure some quenching s ign atures we hope to obs erve. 

(Figure [l4|) and AF (Figure 


Both 




12) cor¬ 


relate with quiescent fraction, while AM* does not. Ihus, 
compactness so defined can be thought of as tracing well 
the average current SFR/M*, but perhaps less well the in¬ 
tegrated effects of quenching on M* over time. Another way 
to say this is that a galaxy’s M*-independent structure may 
provide additional insight into the long-term effects of feed¬ 
back processes in massive galaxies. 


6 DISCUSSION 


Analyzing synthetic data from simulations alongside galaxy 
surveys promises to be a valuable approach to constrain the 
physics of galaxy evolution. Powerful pipelines have been 
developed to convert theoretical models of galaxy forma¬ 
tion into useful mock surveys, such as the Millennium Run 


Database and Observatory (Lemson 2006 Overzier et al. 


2013), the Theoretical Astrophysical Observatory (Bernyk 


et al. 2014), and the simple semi-empirical approach by 
Taghizadeh-Popp et al. (2015). Moreover, such models now 


include treatments motivated by high-resolution simulations 
for the morphological transformations that may be a key 
aspect of galaxy evolution (e.g., Porter et al. 2014[ Bren¬ 
nan et aL]|2015 ). Now, owing to advances in our ability to 
evolve hydrodynamics coupled to sub-grid galaxy formation 
physics, multiple groups have simulated ever-larger volumes 
of a ACDM universe with high enough spatial resolution to 
predict the internal structures of galaxies to varying degrees 
(e.g., Munshi et al.||2Q13 Ceverino et al.|[2014 Christensen 


et al. |2014| [Dubois et al.||2014| |Cen||2014| [Khandai et al!] 


2015| ^gelsberger et al.||2014b| Schaye et ai.|201^ 

An obvious way forward is to begin to combine syn¬ 
thetic observatory pipelines with cosmological hydrodynam- 
ical simulations, because such calculations follow directly the 
gravitational and physical processes which ultimately con¬ 
trol galaxy morphologies. By constraining our models with 
observed galaxy populations, we can make progress toward 
improving not only our understanding of galaxy physics, but 
also our ability to link these processes to observations. 

This requires two actions: creating realistic synthetic 
data and connecting observational diagnostics measured 
from such data to these physical processes. For example, 
this technique has been used to measure the galaxy merger 
rate at z < 1 by measuring how morphologies evolve in 

|2010a|b |. 


merger simulations (Lotz et al. 


2008a 


studies have used cosmological simulations to determine the 
impact of powerful feedback on clumpy high-redshift galax¬ 


ies (Mandelker et al. 2014 Moody et al. 2014a), and how 


galaxy structures might evolve (e.g., Scannapieco et al. 2010 


Pedrosa et al.|20f4 Wellons et al.|2015 ), including by direct 

analogy with HST surveys of distant galaxies (Snyder et al. 


2015). 


To begin expanding this effort to very large-scale hy- 
drodynamical simulations, T15 converted the z = 0 output 
of the Illustris Simulation ( Vogelsberger et al.||2014a Genel 
et al.|2014 ) into synthetic optical images and spectra, a pub¬ 
lic dataset with which it is possible to address a number of 
open scientific questions. Moreover, this observationally ori¬ 
ented approach is an effective way to compare multiple large- 
scale yet moderate-resolution simulation projects, for exam¬ 
ple the Eagle simulation (Schaye et al.||2014 Crain et al 
^OT^, Horizon-AGN ([Dubois et al.||2Q14| ), MassiveBlack- 
II ( Khandai et al.[[2015 ), and others, as well as with future 
larger and higher-resolution calculations. 

In the present work, we measured a suite of non- 
parametric morphology statistics of galaxies in the Illustris 
Simulation at low redshift, and began to explore the extent 
to which these resemble real observations. We have made 
several simplifying choices, including to separate galaxies 
from their neighbors (Section [2.1[ ) and to neglect for now the 
effect of dust on optical morphology. Of course, the quan¬ 
titative morphology of individual galaxies does depend on 
dust ( Lotz et al.|20Q8a Snyder et al.|2015 ), but primarily in 
cases such as starbursts which are rare in the local galaxy 
population. Moreover, we have tested that the optical mor¬ 
phology trends we analyzed are, to first order, unchanged 
by including a treatment for dust. 


6.1 General optical morphologies 

In Section we presented the optical morphology distribu¬ 
tions of Illustris galaxies at = 0. We found that simulated 
galaxies occupy roughly the observed G-M 20 space, with cor¬ 
rect correlations of morphology with stellar mass and SFR. 
This is a noteworthy success: hydrodynamical simulations 
can now produce galaxy populations with reasonable mor¬ 
phology distributions. 

That morphologies occupy the space of SFR, M*, 
galaxy density, and R 1/2 (Figures and 10) each in a re¬ 
alistic fashion was somewhat surprising. The parameters of 
the Illustris physics models were tuned to match observed 
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distributions in both M, and SFR (the axes of Figure 
panel a), but not to recover the median morphology as a 
function of these quantities. Moreover, [Snyder et al. ( 2015| 
showed the dependence of G-M 20 on M* and SFR to be 
moderately insensitive to the precise physics of feedback, 
using very different simulations. 

One explanation is that in a simulation such as Illus- 
tris, a galaxy’s morphology and star formation rate are not 
independent: the rate and mode of the feedback-regulated 
gas supply ( Nelson et al.||2015a ) determines the subsequent 
conversion into stars. If those two things are modeled reason¬ 
ably accurately, then realistic galaxy structures are a natural 
consequence. 

However, another interpretation is to say that the dis¬ 
tribution of galaxy morphology - the coarse breakdown of 
bulges and discs - does not add significantly more infor¬ 
mation to galaxy formation models beyond that provided 
by constraints on the M* and SFR distributions. Indeed, it 
seems that the success of Illustris with the broad features 
of G-M 20 space could be a direct consequence of tuning the 
feedback models to reproduce the stellar mass function and 
stellar halo occupation function at z = 0 and total SFR 
versus time. 

If true, then the regulation of star formation in a crudely 
correct manner, combined with jointly modeling baryons 
and dark matter to accurately follow collapse, cooling, and 
angular momentum transport, would be sufficient to repro¬ 
duce the distribution of automated galaxy morphologies. 
This hypothesis can be tested by implementing more re¬ 
alistic treatments for galaxy physics, by analyzing more de¬ 
tailed morphological and kinematic tracers, and by studying 
directly galaxy structures at fixed DM halo mass or other 


a function of environment in Illustris, similar to the results 


by Cappellari et al. (2011b) 


properties (Section 5.3), where possible. Even if we reach a 
point where these do not change drastically with further im¬ 
provements to galaxy physics, it is clearly desirable to learn 
how the microphysics of star formation and feedback (e.g., 
Hopkins et al.|2014 ) translate into the assumptions used by 
large-scale models (e.g., Muratov et al.||2015 ). 


6.2 Rotation versus optical structure 

In Section 1^ we showed how optical morphology statistics 
relate to the physical kinematics of massive galaxies in Il¬ 
lustris. Broadly, we measure bulge-dominated light profiles 
in galaxies lacking significant rotation, and a diversity of 
light profiles in rotating galaxies with M* > lO^^M©. Al¬ 
though quiescence is correlated tightly with bulge strength 
measured by E(G, M 20 ) (Equation]^, it correlates less with 
kinematic measures of bulge strength and lack of rotation. 
One interpretation of this is that it simply follows recent 
thinking about galaxy structures trigg ered by kinematic sur- 
veys like Sauron and Atlas^^ (e.g., Emsellem et al.||20Q7 


Cappellari et al.|2011a ): flat, rotating galaxies are closely re¬ 


lated, and differ primarily by having differing levels of star 
formation and therefore appear to have stronger or weaker 


spiral arms (e.g., as proposed and discussed by van den 
Bergh|1976 Bender et al.|1994 Kormendy Bender |1 996 ). 
Slowly rotating galaxies are rare, quenched, and occur at 
a similar rate relative to the total quenched fraction inde¬ 
pendent of environment (Section |4.3[ ), as they do in e.g., 
Atlas^^. Also, the relative populations of late-type galax¬ 
ies and fast-rotating early type galaxies evolve smoothly as 


As in real nearby galaxies, Illustris galaxies appear to 
have structures implying that meaningful intrinsic classifi¬ 
cation is more like a “comb” than a “tuning fork” (e.g., 
Eigure 1 by Cappellari et al.||2011b ). Therefore, large-scale 
galaxy formation simulations like Illustris and Eagle are 
capable of producing not only a roughly correct distribution 
of bulge- and disc-dominated galaxies, but also a roughly 
correct phenomenology of galaxy morphology and kinemat¬ 
ics at 2 : = 0. Moreover, Genel et al. (2015) found that the 


angular momentum content of Illustris galaxies depends on 
mass in roughly the same manner as observations and pre¬ 
dictions of analytic models (e.g., Romanowsky Eall|2012 ) 


Therefore, observational constraints which have been tested 


et al. 

2014b 

) or small samples of cosmological simulations 

(Naab et al. 

2014) can also now be tested with large-scale 


hydrodynamical simulations. This will allow us to exploit 
not only the rough diversity of galaxy types, but also their 
detailed statistical distributions, in order to better constrain 
models of galaxy formation. 


6.3 Detailed morphologies 

Although we are finding that simulations can reproduce av¬ 
erage kpc-scale morphologies, we are also beginning to find 
hints of new constraints available in higher-order morphol¬ 
ogy measurements and residual correlations. 

Eor instance, a significant number of Illustris galaxies 
appear to have a ring-like morphology. One possibility is that 
the ring-like shapes result from the equation of state ISM 


pressurization in the Springel & Hernquist (2003) model. 


Another possibility is that the Illustris ISM and feedback 
model is leading to the formation of ring-like star formation 
patterns driven by the choice of density scale at which to ei¬ 
ther form stars or recouple the wind material to the galaxy 
material. In large galaxies this wind mass can recouple at 
a location internal to the galaxy, enhancing the ISM den¬ 
sity and causing additional star formation there (a positive 
feedback effect similar to, e.g., Dugan et al.||2014 ). 

The fact that we can detect this ring-like feature in 
a statistical sample of simulated galaxies, both in young 
and old stars, implies that we may now be able to use it 
as a constraint for large-scale models of galaxy formation. 
The number densities, masses, ages, SERs, and sizes of ring¬ 
like galaxies may rule out or better constrain the numeri¬ 
cal treatment of the ISM and feedback models like the one 
used in Illustris. This could require treating dust faithfully 
in these images. 


6.4 Morphologies at intermediate and low mass 

There is a lack of quenching in Illustris galaxies with M* < 
IqIo.s^©, both to stellar mass functions that are too 


high and also quenched fractions that are too low (Vogels- 
berger et al. [201481 Genel et al.|2Q14 ). This appears also as 
an apparent shortfall in the numbers of quenched spheroids 
in the distribution of morphology versus M* and SER at 
intermediate masses, such as we see at M* ~ 
in Eigure 1^ (Sections 3.2 and 3.3). One obvious interpre- 
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tat ion is that feedback is too weak in these systems. How¬ 
ever, the morphologies in low mass galaxies may also be 
affected by other factors, such as the inability to follow gas 
ffows toward galaxy centers on very small scales, a mech¬ 
anism which could cause the remnants to contract. Also, 
the artificial ISM pressurization used by Illustris and other 
simulations of comparable volume could prevent gas from 
being influenced by turbulence and collapse on very small 
physical scales <C 1 kpc. This could prevent fragmentation 
and inward migration of gas during periods of “violent disc 
instability” (VDI), which could assist the formation of com¬ 
pact spheroids (e.g., |Dekel et al.|2009 Zolotov et al.|2015 ) in 
intermediate mass galaxies. Also, as in Illustris, lower-mass 
halos (Mh < 10 ^^Mq) in both semi-analytic and large-scale 
hydrodynamical models overproduce stars at high redshift 
( White et aL][M^ . 


6.5 Residual morphology correlations 


If multiple galaxy formation models can broadly reproduce 
the diversity of galaxy types, then we should appeal to resid¬ 
ual or higher-order correlations in order to make new falsifi- 
able predictions. The Illustris galaxy physics model includes 
a component based on feedback from SMBHs, which im¬ 
prints a residual dependence on morphology at fixed 


and Ri /2 (Section 5.2). Sijacki et al. (2015) presented the 
prediction for this residual correlation with galaxy colours. 
Similarly, this model imprints residual correlations in 


galaxy morphologies at fixed halo mass (see also Pillepich 
|et al.|20l4 ). In Sectionwe found that at fixed halo mass 
(for 12 < log^QM/i/M© < 12.5), M* is inversely propor¬ 
tional to optical bulge strength as measured by T(G, M 20 ) 
(Equation]^. In other words, in galaxies at fixed halo mass, 
more disc-dominated systems have greater stellar mass than 
more bulge-dominated ones. Relatedly, at fixed halo mass, 
M* is inversely proportional to Mbh, so that in galaxies 
with fixed halo mass, over-massive SMBHs imply less total 
stellar mass than in galaxies with under-massive SMBHs. 
Necessarily then, Mbh and T(G, M 20 ) are positively corre¬ 
lated at fixed Mh- 

This appears to be a natural consequence of the physics 
of feedback assumed in Illustris. In a given halo, more SMBH 
accretion ultimately deposits more energy into the galaxy’s 
gas reservoir, preventing it from forming as many stars as 
it otherwise would. Since stars form most often in cold gas 
that has collapsed into a disc, this quenching signature is 
also imprinted in galaxy morphology: more quenching leads 
to fewer young stars and a fainter disc component, leading 
eventually to a lower M* and earlier type in galaxies of the 
same halo mass. This reduction in M* at fixed Mh may 
imply that we should contrast the behavior of M*-based 
structural tracers (e.g., some measures of compactness) with 
M*-independent ones like G-M 20 , Sersic index, and C, for 
example. 

This hypothesis may be difficult to disentangle from al¬ 
ternative ones. For example, what causes a SMBH to be 
over-massive for a given Mh in the first place? A common 
answer might invoke processes leading to the co-evolution of 
bulges and SMBHs, in which case galaxy morphology vari¬ 
ations might precede quenching variations at fixed Mh- To 
find out what Illustris predicts for this question, it will be 
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important to map how individual simulated galaxies evolve 
over time in distributions of these important properties. 


7 SUMMARY AND CONCLUSIONS 

We converted z = 0 galaxies formed in the Illustris simula¬ 
tion into synthetic images and measured their optical mor¬ 
phology using automated tools. From idealized images, we 
created 42531 SDSS-like sources for 10808 simulated galax¬ 
ies (up to 4 angles per galaxy) with lO^'^M© < M* < 
lO^^'^M©, in each of rest-frame u, g, z, and H filters. From 
these synthetic data, we measured non-parametric diagnos¬ 
tics of galaxy morphology, catalogs of which we present in 
Appendix [b| In equation we defined a non-parametric 
bulge-strength parameter, F{G, M 20 ), which is easy to mea¬ 
sure and relatively robust against the effects of dust and 
disturbances. 

Our main conclusions about the morphology of galaxies 
formed by the Illustris Simulation at z = 0 are: 

(i) The distribution of synthetically observed morphol¬ 
ogy agrees well with the observed one at z ~ 0. Crucially, 
this includes a substantial population of quenched bulge- 
dominated galaxies having C ~ 0.6 and M 20 ~ —2.5, and 
a locus of disc-dominated galaxies with realistic concentra¬ 
tions and pixel light distributions. 

(ii) At M* ~ lO^^M©, we find that Illustris galaxy struc¬ 
tures are a combination of spirals and composite systems 
with G < 0.5 and M 20 > —2. As M* increases above 
10^^M©, the primary locus shifts smoothly toward the re¬ 
gion occupied by ellipticals. 

(hi) At 10^°'^ < M^IMq < 10^^, simulated G-M 20 values 
are bimodal, with a separate peak at G ~ 0.4 and M 20 > — 1. 
The extreme location of this second peak is a consequence of 
a large number of galaxies having a “ring-like” morphology. 
A detailed comparison of the population of such galaxies in 
Illustris with observed samples may be a new constraint on 
models of large-scale galaxy formation. 

(iv) As a function of M* and SFR, the simulated optical 
morphologies follow trends observed in galaxy surveys: at 
fixed M*, lower-SFR galaxies have earlier structural types, 
and at fixed SFR, higher-mass galaxies have earlier struc¬ 
tural types. 

(v) Optically defined morphologies correlate with kine¬ 
matically defined ones, especially in simulated galaxies with 
M* > 10^^M©. Some kinematically defined tracers are less 
correlated with quiescence than is G-M 20 . 

(vi) Simulated optical morphologies trace the same 
rot at ion-SFR-environment relation as found by integral field 
spectroscopy surveys such as Atlas^^: slowly rotating, 
quenched ellipticals are rare and occur as a constant share of 
total quenched galaxies. Rotating galaxies are diverse, some 
with high SFR and spiral arms while others are “anemic 
spirals”: flat, rotating, but quenched. 

(vii) At M* ~ 10^^ M©, the ^-band size distribution of 
Illustris Simulation galaxies roughly matches the observed 
one. At M* ~ 10^°M©, the Illustris sizes are a factor of ~ 2 
larger than observed. 

(viii) At fixed halo mass Mh and at fixed relative rotation, 
quiescence correlates strongly with optical bulge strength. 

(ix) At fixed intermediate halo mass {Mh ~ 10^^M©), M* 
correlates inversely with optical bulge strength F(G, M 20 ) 
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and SMBH mass, which we believe results from the Illustris 
AGN feedback model. In other words, AGN quenching acts 
to reduce the total stellar mass below what would have oth¬ 
erwise formed, with a side effect of reducing the extent to 
which the light profile is disc-dominated. 

The Illustris simulation is a large-scale calculation of 
galaxy formation with physics model parameters set to ap¬ 
proximately reproduce selected global statistics. These in¬ 
clude the galaxy stellar mass function, the stellar mass-halo 
mass relation, and the evolution of total star formation den¬ 
sity. The field has made significant efforts to overcome nu¬ 
merical challenges that had previously limited the ability 
of cosmological simulations to form and maintain a diverse 
population of star-forming disks and quiescent spheroids. 

We have shown that these ingredients lead to a di¬ 
verse distribution of galaxy structures that includes both 
spheroids and disks at z = 0. The non-parametric mor¬ 
phologies of our simulated galaxies agree well with observed 
populations. Gertainly the balance of these populations is 
imperfect in Illustris, since, for example, quenching is not as 
complete as observed in galaxies with M* ~ lO^^M©. This 
leads in part to the simulated galaxies having larger radii, 
a less-clear structural bimodality, and more extended disks 
than observed. However, we find an approximately correct 
correlation between star formation and galaxy morphology, 
which we believe merits optimism. It implies that the ingre¬ 
dients above, plus physics improvements and/or higher spa¬ 
tial resolution, could be sufficient to reproduce a more fully 
realistic galaxy population. In other words, we hypothesize 
that for a model which properly regulates a galaxy’s star 
formation, and which mitigates possible impacts on second- 


order structural parameters (see Section 3.3), structure for¬ 


mation plus galaxy physics leads naturally to the distribu¬ 
tion of galaxy morphologies. Proving this hypothesis will 
require comparing the several recent large-volume hydrody- 
namical simulations, performing a parameter study of large- 
volume simulations with different feedback implementations, 
or both. 

We have also presented our methods for modeling syn¬ 
thetic observations of galaxies. This includes the release of 
Python modules for the manipulation of idealized synthetic 
images (in our case, output of the Sunrise code) into ob¬ 
servational analogues, which can then be analyzed with au¬ 
tomated tools and by visual classifiers. 

We expect to use these methods to enable additional 
science with galaxy surveys. With the emergence of statisti¬ 
cally relevant cosmological simulations that directly include 
much of the required galaxy physics, we can, for instance, 
measure the predicted observability time of key galaxy for¬ 
mation phenomena, and improve the robustness of the nec¬ 
essary galaxy classification schemes. Moreover, T15 demon¬ 
strated several projects based on detailed simulated images, 
including how shells observed in massive spheroids might re¬ 
flect their formation histories, how SED fitting routines can 
be tested against simulated SED maps, and generally how 
the spectra of the galaxy population could be used to fur¬ 
ther constrain galaxy formation physics. Such calculations 
will improve measurements of the galaxy-galaxy merger rate 
from disturbed morphologies, the formation of disks and 
bulges, the impact of feedback, and the role of star-forming 
clumps. 
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APPENDIX A: LIGHT ASSIGNMENT TESTS 

Here we show how our method of assigning sizes to star par¬ 
ticles affects the morphology measurements studied in this 
paper. In Sunrise, star particles emit light with a distri¬ 
bution approximating a simple SPH kernel with smoothing 


length r* (Jonsson 2006). Our fiducial method assigns r* 
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equal to the radius which encloses a star particle’s 16 near¬ 
est star particle neighbors (T15), which we denote rl^. The 


most 

common alternative method (e.g., Scannapieco et al. 

2010 

Bush et al.|2010| Snyder et al.|2011) assigns star par- 


tides a constant radius equal to the gravitational softening 
length for stars; in Illustris at z = 0, r* = 0.710 kpc. 


Figure |A 1 | compares these two methods as a function 
of PSF FWHM for five example galaxies at various masses 
(9.5 < M*/M 0 < 11.5) placed at z = 0.05 and a pixel size of 
0.25 arcsec. Since we model observations with PSF FWHM 
near the Illustris-1 simulation spatial resolution, which sets 
the star particle sizes in either method, these images are 
not especially sensitive to PSF size. We hnd that using a 2 
arcsec PSF FWHM instead of 1 arcsec changes T(G, M 20 ) 
(Equation the quantity used throughout this paper) by 
a median (mean) of —0.01 (0.03), with a median absolute 
difference of 0.06. 

In this paragraph, we fix PSF FWHM= 1.0 arcsec ~ 1 
kpc, the value we employed in this paper. Across all five 
galaxies with four viewing angles each, the median (mean) 
difference in M 20 between the two light assignment methods 
is 0.009 (0.026), or 0.85% (0.89%). The median (mean) dif¬ 
ference in G is 0.006 (0.005), or 1.0% (1.5%). The median 
(mean) difference in F{G^ M 20 ) is 0.03 (0.01), with a median 
absolute difference of 0.06. The sense of this difference is that 
on average, image pixels are slightly less concentrated (M 20 
is larger) but more unequal {G is larger) under our adap¬ 
tive assumption versus the constant r'^ assumption. The 
magnitude is small enough to have no effect on our conclu¬ 
sions. 


APPENDIX B: MORPHOLOGY CATALOGS 

We make the morphology measurements from this paper 
available as supplementary tables associated with this arti¬ 
cle. These catalogs will also be publicly available on the Illus¬ 
tris Project website, illustris-project.org. Moreover, the code 
we used to add observational effects to the images of T15 
is publicly available at https://bitbucket.org/ptorrey/sunpy. 
The ideal synthetic images are already publicly available 
(T15). 

Table [BT] presents example entries of these catalogs. We 
imaged each of the 10808 simulated galaxies with M* > 
lO^'^M© from four viewing directions in four filters. Since 
these all lead to different quantitative morphologies, each of 
these 16 iterations per simulated galaxy is treated as a sepa¬ 
rate source or line in our catalogs. We separate our tables by 
hlter, leading to 4 individual catalogs of roughly 43000 en¬ 
tries. Filesystem errors cause a very small number (< 1%) of 
the morphology measurements to fail at random, and so we 
have cross-matched our catalogs to a uniform set of 42531 
sources per hlter. Finally, we sort by M* and show the top 
20 lines of our unattenuated ^-band catalog in Table |B1| 
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Figure Al. Face-on example ^f-band images from Illustris at 2 ; = 0, as a function of star particle radius assignment method and 
PSF FWHM, labeled with G and M 20 values at the top left and top right of each panel, respectively. The blue solid bordered column 
represents the images used throughout this paper, while the orange dotted border shows the most common alternative (fixed r'^ = 0.71 
kpc). 


Table Bl. First 20 entries of our unattenuated gf-band morphology catalog, sorted by M*. The full catalogs for n, z, and H will be 
made available as online-only supplementary material. The subhalo index is the index into the Subfind FOF catalog arrays for the 2 ; = 0 
snapshot (number 135). U — B \s the intrinsic unattenuated AB-system colour of the galaxy’s stellar populations. The camera index is 
the Sunrise camera number. The remaining quantities are morphology measurements described in Section |2.3| 
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